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We present results on the chiral and deconfinement properties of the QCD transition at finite 
temperature. Calculations are performed with 2 + 1 flavors of quarks using the p4, asqtad and 
HISQ/tree actions. Lattices with temporal extent Nt = 6, 8 and 12 are used to understand and 
control discretization errors and to reliably extrapolate estimates obtained at finite lattice spacings 
to the continuum limit. The chiral transition temperature is defined in terms of the phase transition 
in a theory with two massless flavors and analyzed using 0{N) scaling fits to the chiral condensate 
and susceptibility. We find consistent estimates from the HISQ/tree and asqtad actions and our 
main resuh is Tc = 154 ± 9 MeV. 
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I. INTRODUCTION 



It was noted even before the advent of Quantum Chromodynamics (QCD) as the underlying theory of strongly 
interacting elementary particles, that nuclear matter cannot exist as hadrons at an arbitrarily high temperature or 
density. The existence of a limiting temperature was formulated in the context of the Hagedorn resonance gas model 
This phenomenon has been interpreted in the framework of QCD as a phase transition ^ separating ordinary 
hadronic matter from a new phase of strongly interacting matter — the quark gluon plasma [3| . Today, the structure 
of the QCD phase diagram and the transition temperature in the presence of two light and a heavier strange quark 
is being investigated using high precision simulations of lattice QCD. 

Understanding the properties of strongly interacting matter at high temperatures has been a central goal of nu- 
merical simulations of lattice QCD ever since the first investigations of the phase transition and the equation of state 
in a purely gluonic SU(2) gauge theory Early work showed that chiral symmetry and its spontaneous break- 

ing at low temperatures play an important role in understanding the phase diagram of strongly interacting matter. 
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Chiral symmetry breaking introduces a length scale, and the possibility that it may be independent of deconfinement 
phenomena was discussed 0, [3- Similarly, the consequences of the existence of an exact global symmetry in the 
chiral limit of QCD, the spontaneous breaking of this 0(4) symmetry, the influence of the explicit breaking of the 
axial ?7a(1) symmetry and the presence of a heavier strange quark for the QCD phase diagram were analyzed 0]. 
The possibilities that the strange quark mass could be light enough to play a significant role in the QCD transition, 
and/or an effective restoration of axial symmetry may trigger a first order phase transition in QCD for even nonzero 
light quark masses were also discussed. Neither of these situations seems to be realized in QCD for physical light and 
strange quark masses. Based on recent high precision calculations, the QCD transition at nonzero temperature and 
vanishing chemical potentials is observed to be an analytic crossover |lC)| - [l3 |. 

In this paper we focus on understanding the universal properties of the QCD phase transition in the chiral limit and 
extract the behavior of QCD at physical quark masses using an 0{N) scaling analysis. We also calculate quantities 
that probe the deconfinement aspects of the QCD transition: quark number susceptibilities and the renormalized 
Polyakov loop. It is expected that for temperatures below the transition temperature these quantities should be well 
described by the hadron resonance gas (HRG) model which is very successful in describing thermodynamics at low 
temperature and the basic features of the matter produced in relativistic heavy ion collisions [isj . It is, therefore, 
interesting to quantify the interplay between the universal properties of the chiral transition and the physics of the 
HRG model. Some of these issues will be addressed in a separate publication. 

Investigations of QCD at finite temperature are carried out using a number of different lattice formulations of the 
Dirac action. While studies based on the Wilson [l3.[l5j or chiral fermion formulations are, at present, constrained 
to a regime of moderately light quark masses (m;/ms>0.2), calculations exploiting staggered fermion discretization 
schemes [lol [Til [T7l - [2^ can be performed with an almost realistic spectrum of dynamical light and strange quarks. 
Today, high statistics calculations, performed at a number of values of the lattice cutoff and quark masses, allow for 
a detailed analysis of discretization errors and quark mass effects. 

Recent studies of QCD thermodynamics with two degenerate light quarks and the heavier strange quark have been 
performed with several staggered fermion actions that differ in the way improvements are incorporated to reduce 
the effects of known sources of discretization errors. These include the asqtad, p4 and stout actions. The results at 
a 0.1 fm show differences not only in the determination of relevant scales, such as the transition temperature, but 
also in the temperature dependence of thermodynamic observables. Because estimates of all observables ought to 
agree in the continuum limit, discretization errors and the dependence on light quark masses require careful analysis. 
Furthermore, since QCD for physical quark masses docs not display a genuine phase transition, the definition of the 
transition temperature itself requires care. A proper definition of a pseudocritical temperature should be related to 
the chiral phase transition in the massless limit of QCD and reduce to it in that limit. The aim of this paper is to 
study chiral and deconfinement aspects of the QCD transition at sufficiently small lattice spacing as to give control 
over the continuum extrapolation, and demonstrate the consistency of the results obtained with different actions. 
Furthermore, for the first time we provide a determination of the chiral transition temperature in the continuum limit 
that makes close connection with the critical behavior of QCD for massless light quarks. 

This paper extends earlier calculations, performed with the asqtad and p4 actions on lattices with temporal extent 
Nr ~ 8 and light to strange quark mass ratio mi/uis = 0.1, in several ways. We have added calculations for 
mi /ma = 0.05 on Nt ~ 8 lattices and performed new calculations with the asqtad action at smaller lattice spacing, 
i.e., for N-r = 12 with mi/ms = 0.05. More importantly, we have performed thermodynamic calculations with the 
highly improved staggered quark action (HISQ) t25i] on lattices with temporal extent A^,- = 6, 8 and 12 to quantify 
discretization errors. Preliminary versions of the results given in this paper have been presented in Refs. |26l-l34l|. 

This paper is organized as follows. In the next section, we discuss thermodynamics calculations with improved 
staggered fermion actions, in particular, emphasizing the HISQ formulation which has been exploited in this context 
for the first time. We analyze the so-called taste symmetry violations in different improved staggered actions. Details 
of the simulation parameters used in our calculations are also given in this section. In Sec. IIIIl we introduce the basic 
observables used in the analysis of the QCD transition and discuss their sensitivity to the expected critical behavior in 
the chiral limit. We present our numerical results for chiral observables in Sec. IIVI In Sec. El we discuss the universal 
properties of the chiral transition and the determination of the pseudocritical temperature. The deconfining aspect 
of the QCD transition, which is reflected in the temperature dependence of the quark number susceptibilities and 
the renormalized Polyakov loop, is discussed in Sec. IVII Finally, Sec. IVIII contains our conclusions. Details of the 
simulations, data and analysis are given in the appendices. 
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II. THERMODYNAMICS WITH STAGGERED FERMION ACTIONS 
A. Staggered fermion actions 

All staggered discretization schemes suffer from the well known fermion doubling problem, i.e., a single staggered 
field describes four copies of Dirac quarks. These extra degrees of freedom are called taste and the full taste symmetry 
(degeneracy of the four tastes) is realized only in the continuum limit. At nonzero lattice spacing a taste symmetry is 
broken and only a taste non-singlet axial U{1) symmetry survives. Consequently, in the chiral limit there is a single 
Goldstone meson, and the other 15 pseudoscalar mesons have masses of order a^a^. For lattice spacings accessible in 
current numerical studies, the effects of taste breaking can produce significant distortions of the hadron spectrum. In 
thermodynamic calculations, these effects are expected to be most significant at low temperatures where the equation 
of state is governed by the spectrum of hadrons (ssj. Current lattice calculations show that the distortion of the 
spectrum accounts for a large part of the deviations of the QCD equation of state from a hadron resonance gas 
estimate [35l - [38| . 

Several improvements to the staggered fermion formulation have been proposed to reduce the O(a^ ) ta ste symmetry 
breaking effects. These improvements involve using smeared gauge fields or the so-called fat links [39[ by including 
paths up to length seven in directions orthogonal to the link being fattened. With these improvements it is possible 
to completely cancel taste symmetry breaking effects at order aga'^ 0. Such an action is called asqtad and has been 
studied extensively [4l|. Fat links are the sum of SU{3) matrices corresponding to different paths on the hyper-cubic 
lattice, and are not elements of the SU{3) group. It has been shown that projecting the fat links back to SU{3) 
(4^ [43 or even to the U (3) grou p [4^ greatlv improves the taste symmetry. Projected fat links are being used in 
simulations with the stout action |22l424l | and HISQ action [2^, |4^, |46| . In this paper, we confirm that this projection 
results in reductions of taste symmetry violations and a much better reproduction of the physical hadron spectrum 
in calculations starting at moderately coarse lattice spacing, i.e., a ~ 0.15 fm. 

In the generation of background gauge configurations, the reduction of the number of tastes from four to one for 
each flavor uses the so-called rooting procedure, i.e., the fermion determinant in the QCD path integral is replaced 
by its fourth root. Effectively this amounts to averaging over the non-degenerate spectra of mesons and baryons, e.g., 
over the non-degenerate spectrum of sixteen taste pions. The validity of this procedure is still a subject of debate 
(47I [isj . (For a more detailed summary of the issues, see Ref. [llj.) Reducing taste symmetry violations is, in any 
case, important for making the rooted staggered theory a good approximation to a single flavor physical theory. 

In the context of thermodynamic calculations, it is also important to control cutoff effects that manifest themselves 
as distortions of the high temperature ideal gas and perturbative high temperature limits. To reduce these 0{a^) effects 
we use improved staggered fermion actions that include three-link terms in the discretization of partial derivatives in 
the Dirac action. These three-link terms remove the tree-level 0{a^) discretization effects, which arc the dominant 
ones at high temperatures j49l[50| . as can be seen by considering the free energy density in the ideal gas limit calculated 
on four dimensional lattices with varying temporal extent Nt. This free energy density of a quark gas divided by the 
corresponding result in the continuum limit (Nr — > 00) is shown in Fig. [TJ For the unimproved staggered fermion 
action with 1-link discretization, there is a significant cutoff dependence for N^- < 16. Including three-link terms in 
the action (p4 and Naik) reduces cutoff effects to a few percent even for Nt = 8. The Naik action with straight 
three-link terms is the building block for both the asqtad and HISQ actions; however, projection of fat links to f/(3) 
to further reduce taste violations is done only in the HISQ action. The stout action, on the other hand, uses just the 
standard 1-link discretization scheme with stout smeared links that include projection to SU{3). 

The HISQ action improves both taste symmetry breaking [2^ and cutoff effects in the hadron spectrum which, as 
mentioned above, are of particular relevance to thermodynamic calculations at low temperatures. The construction 
of the projected fat link action proceeds in three steps. In the first step, a fat7 link is constructed; i.e., a fat link 
which includes all the paths in orthogonal directions up to length seven. This step is common to the asqtad action. 
In step two, the sum of the product of SU{3) matrices along these paths is projected to U{3). In the third step, these 
projected fat links are used in the conventional asqtad Dirac operator without tadpole improvement. Thus, from the 
point of view of reducing taste symmetry breaking at order 0{asa^) the asqtad and the HISQ actions are equivalent, 
but differ at higher orders. Unfortunately, these higher order terms are large in the asqtad action as discussed in 
Sec. HID I where we show that the projection of fat links to U{3) in the HISQ formulation significantly reduces the 
distortion of the spectrum at low temperatures. The straight three-link Naik term in the asqtad and HISQ actions 
eliminates the tree-level C(a^) discretization effects, consequently, their behavior at high temperatures is equivalent. 

For the HISQ calculations presented here, we use a tree-level improved Symanzik gauge action that is also common 
to the p4 and stout formulations. We refer to this combination of the gauge and 2-t-l HISQ quark actions as the 
HISQ/tree action to distinguish it from the HISQ action used by the MILC collaboration in their large scale zero- 
temperature 24-14-1 fiavor simulations with a dynamical charm quark [45|,|4^. In the 2-t-l-|-l HISQ action, in addition 
to the 1-loop tadpole improved version of the Symanzik gauge action, the 1-loop and mass dependent corrections are 
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Figure 1: The free energy density of an ideal quark gas calculated for different values of the temporal extent Nr divided by the 
corresponding value for Nt = oo. 



included in the Naik term for the charm quark. 



B. Lattice Parameters and Simulation Details 



A summary of the run parameters, statistics, and data for the p4, asqtad, and HISQ/tree actions analyzed in this 
paper is given in Appendix [A] We have previously presented the equation of state and other thermodynamic quantities 
using the p4 and asqtad actions on lattices vifith temporal extent Nr = 4, 6, and 8 in Refs. [l8l - [2l| . Here, we extend 
these studies in the following three ways: 

• Additional asqtad calculations for N^- = 8 with mi /rag = 0.2 and 0.05. (See Table Hill in Appendix [AD 

• New simulations with the asqtad action in the transition region on lattices with temporal extent Nt = 12 and 
light to strange quark mass ratio mi/ms = 0.05. (See Table IV] in Appendix H) 

• New results with the HISQ/tree action on N^. = 6 lattices with mj/m^ = 0.2, 0.05 and 0.025; on N^ = 8 lattices 
with mi/ms = 0.05 and 0.025; and on Nr = 12 lattices with mz/m^ = 0.05. (See Tables IVl IVTl and IVlIl in 
Appendix E)) . 

As described in previous studies, the first step is to determine the line of constant physics (LCP) by fixing the strange 
quark mass to its physical value mg at each value of the gauge coupling /? (20l.[2l|^. In practice, we tune rus until the 
mass of the fictitious rjss meson matches the lowest order chiral perturbation theory estimate m,,^5 = •\/2m'^ — 
(lol Hlj. Having fixed m^, between one and three values of the light quark mass, mi/ms = 0.2, 0.1, 0.05 and 0.025 
are investigated at each Nr and used to obtain estimates at the physical point mi/mg = 0.037 by a scaling analysis 
discussed in Sections IIIII and HVl 

All simulations use the rational hybrid Monte-Carlo (RHMC) algorithm [5l|,[52] ■ The length of the RHMC trajectory 
is 0.5 in molecular dynamics (MD) time units (TU) for the p4 action and 1.0 for the asqtad action. The statistics, 
therefore, are given in terms of TUs in Tables ITll - IVIII in Appendix |Al For each value of the input parameters, we 
accumulated several thousand TUs for zero-temperature ensembles and over ten thousand TUs for finite temperature 
runs. The RHMC algorithm for the HISQ action is discussed in Refs. [1^, [s^. In the calculations with the HISQ/tree 
action, the length of the RHMC trajectory is typically one TU. For smaller values of (coarse lattices) we used 
trajectories with length of 1/2 and 1/3 TU since frequent spikes in the fermion force term |53| reduce the acceptance 
rate for longer evolution times. 

To control finite size effects, the ratio of spatial to temporal lattice size is fixed at N„/Nr = 4 in most of our finite 
temperature simulations. The exception is Nr = 6 runs at milrris = 0.2, which were done on lattices of size 16"^ x 6. 



In the case of the p4 action and for gauge couplings j3 close to the transition region, the bare quark masses were not varied with /3 to 
facilitate the Ferrcnberg-Swendscn re-weighting procedure. 
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Figure 2: The static potential calculated for the HISQ/tree action with mj = 0.2ms (left) and mi = 0.05ms (right) in units of 
ro- In the plot on the left, we compare the HISQ/tree and p4 results obtained at a similar value of the lattice spacing. The 
dashed line in the plot on the right is the string potential Vstring{r) = —n/{12r) + or matched to the data at r/ro = 1.5. 



Zero-temperature calculations have been performed for different lattice volumes (see Table IVII in Appendix \^ such 
that the spatial extent of the lattice L satisfies LAf^ > 3, except for the smallest lattice spacings where LMt^ ~ 2.6. 

The construction of renormalized finite temperature observables requires performing additive and multiplicative 
renormalizations. We implement these by subtracting corresponding estimates obtained at zcro-tcmperaturc. These 
matching zero-temperature calculations have been performed at several values of the parameters and then fit by 
smooth interpolating functions over the full range of temperatures investigated. 

In the following two subsections, we discuss the determination of the lattice spacing and the LCP for the HISQ/tree 
action, and the effect of taste symmetry breaking on the hadron spectrum. 



C. The static potential and the determination of the lattice spacing 



The lattice spacing is determined using the parameters tq and ri, which are fixed by the slope of the static quark 
anti-quark potential evaluated on zcro-tcmpcrature lattices as [s^ 

\ ^' / r=ro V / r=ri 

and set the scale for all thermodynamic observables discussed in this work. The calculation of the static potential, 
ro and ri for the p4 action was discussed in Refs. [l^ [2l[. In particular, it was noticed that for the values of /3 
relevant for the finite temperature crossover on A^,- = 8 lattices, the parameter rg is the same, within statistical errors, 
for rai = Clm^ and mi = 0.05ms. Therefore, we use the interpolation formula for rg given in Ref. [l9| to set the 
temperature scale for the p4 data. The calculation of the static potential and r\ for the asqtad action was discussed 
in Ref. [4l| . In Appendix [B] we give further details on the determination of ■ Here we note that the statistical 
errors in the rxja determination are about 0.2% for gauge couplings relevant for the AT^ = 8 calculations and about 
0.1% for the A^^ = 12 calculations. We also reevaluate systematic errors in the determination of ri/a and find that 
these errors are smaller than 1% on A^t- — 12 and about 1% on A^,- = 8 lattices. These uncertainties will impact the 
precision with which the chiral transition temperature is estimated. 

The static quark potential for the HISQ/tree action is calculated using the correlation functions of temporal Wilson 
lines of different length evaluated in the Coulomb gauge. The ratio of these correlators, calculated for two different 
lengths, was fit to a constant plus exponential function from which the static potential is extracted. To remove the 
additive UV divergence, we add a /3-dependent constant c(/3) defined by the requirement that the potential has the 
value 0.954/ro at r = rg. This renormalization procedure is equivalent to the normalization of the static potential to 
the string potential Vstringiy) = — 7r/(12r) -l-ar at r = 1.5ro [20|. The renormalized static potential, calculated for the 
HISQ/tree action for mi/ms = 0.05, is shown in Fig. [^^right) and wc find no significant dependence on (3 (cutoff). We 
conclude that discretization errors, including the effects of taste symmetry violations, are much smaller in the static 
potential compared to other hadronic observables. Furthermore, for approximately the same value of rg/a, the static 
potentials calculated with the HISQ/tree and p4 actions agree within the statistical errors as shown in Fig. [2l^left). 

To determine the parameters tq and ri, we fit the potential to a functional form that includes Coulomb, linear, and 
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Figure 3: The ratio ro/ri for the HISQ/tree action. Fitting all the data at /3 > 6.423 by a constant gives ro/ri — 1.508(5) as 
our best estimate of the continuum extrapolated value. 



In this Ansatz, the Coulomb part is corrected for tree-level lattice artifacts by introducing a fourth parameter A. 
The term proportional to A reduces systematic errors in the determination of tq and ri due to the lack of rotational 
symmetry on the lattice at distances comparable to the lattice spacing. The resulting fit has a x^/dof close to unity 
for r/a > -v/S, except for the coarsest lattices corresponding to the transition region on ~ 6 lattices. Consequently, 
we can determine rg/a reliably for /S corresponding to the transition region for N^- = 8 and = 12 lattices. On 
lattices corresponding to the transition region for = 6 we use a three parameter fit (Coulomb, linear and constant) 
with lattice distance replaced by tree- level improved distance, r — > r/. Following Ref. [56j . r/ was determined from 
the Coulomb potential on the lattice. Neither the three-parameter nor the four-parameter fit gives acceptable x^; 
however, the difference in the tq values obtained from these fits is of the order of the statistical errors. We use this 
difference as an estimate of the systematic errors on coarse lattices. 

One can, in principle, extract r^/a and ri/a by using any functional form that fits the data in a limited range about 
these points to calculate the derivatives defined in Eq. ([T]). We use the form given in Eq. but perform separate 
fits for extracting ro/a and ri/a for each ensemble. The fit range about r^/a (or ri/a) is varied keeping the maximum 
number of points that yield x^/dof ~ 1. The variation in the estimates with the fit range is included in the estimate 
of the systematic error. 

The value of ri/a is more sensitive to the lattice artifacts in the potential at short distances than ro/a. Only for 
lattice spacings corresponding to the transition region on Nr = 12 lattices are the lattice artifacts negligible. Again, 
we used the difference between four and three parameter fits to estimate the systematic error in ri/a. 

The lattice artifacts due to the lack of rotational symmetry play a more pronounced role in the determination of 
ri/a, and are significantly larger for the HISQ/tree than for the asqtad action. This is presumably due to the lack of 
tadpole improvement in the gauge part of the HISQ/tree action. For this reason, we use rg on coarser lattices and ri 
on fine lattices and connect the two using the continuum estimate of tq/ti for estimating the scale. Further details of 
this matching are given in Appendix [BJ As noted above, these effects are no longer manifest for the crossover region 
for Nr = 12 lattices. 

Having calculated ro/a, ri/a and r^/ri at a number of values of /?, we estimate the continuum limit value for the 
ratio To/ri. In Fig. |3] we plot the data for the HISQ/tree action. It shows no significant variation with (3, so we 
make two constant fits to study the dependence on the range of points included. The first fit includes all points with 
/3 > 6.423 and the second with /3 > 6.608. We take r^/ri = 1.508(5) from the first fit as our best estimate since it 
has a better 'X^ / dof ~ 0.32, includes more points and matches the estimate from a fit to all points. This estimate is 
higher than the published MILC collaboration estimate using the asqtad action: ro/ri = 1.474(7) (18) [sHi- A more 
recent unpublished analysis, including data at smaller lattice spacings, gives ro/ri = 1.50(1) for the asqtad action [STf . 
consistent with the HISQ/tree estimate. We will, therefore, quote ro/ri = 1.508(5) as our final estimate. 

Finally, to extract the values of ri or tq in physical units, one has to calculate these quantities in units of some 
observable with a precisely determined experimental value. We use the result ri = 0.3106(8)(18)(4) fm obtained 
by the MILC collaboration using to set the lattice spacing [H^. This estimate is in good agreement with, and 
more precise than, the recent values obtained by the HPQCD collaboration: ri = 0.3091(44) fm using bottomonium 



constant terms [ill. Issj: 




(2) 
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splitting, fi = 0.3157(53) fm using the mass splitting of Dg and rjc mesons, and ri = 0.3148(28) (5) fm using fss, the 
decay constant of the fictitious pseudoscalar ss meson [59j . To set the scale using tq, we use the above result for tq/ti 
to convert from ri to rg. This gives rg = 0.468(4) fm, which is consistent with the estimates tq = 0.462(11)(4) fm by 
the MILC collaboration ^ and tq = 0.469(7) fm by the HPQCD collaboration 

There are two reasons why we prefer to use either rg or ri to set the lattice scale. First, these are purely giuonic 
observables and therefore not affected by taste symmetry breaking inherent in hadronic probes. Second, as discussed 
above, ro (and ri) do not show a significant dependence on mi/nis and thus one can extract reliable estimates for 
the physical LCP from simulations at mi /rag = 0.05. Nevertheless, we will also analyze the data using fx to set the 
scale and discuss its extraction in Sec. IIIDl 

D. Hadron masses and taste symmetry violation 

Precision calculations of the hadron spectrum have been carried out with the asqtad action in Refs. [4l|, [5^ [6ll |. 
Details of the calculations of hadron correlators and hadron masses used in this paper are given in Appendix [C] For 
completeness, we also list there the masses of baryons estimated at the same lattice parameters. 

As describe d in Sec. IIIB] the strange quark mass is fixed by setting the mass of the lightest ss pseudoscalar to 
-\/2M|. — M'^ = 686 MeV~[5^. Masses of all other pseudoscalar mesons should then be constant along the lines of 
constant physics defined by mi /nig = 0.2 and 0.05. Fits to the data give 

roM^ = 0.3813(12), tqMk = 1.1956(33), rgM^^, = 1.6488(46), m, = O.OSm^, (3) 
roM^ = 0.7373(14), toMk = 1.2581(23), rgM^^, = 1.6206(30), m/ = 0.20ms. (4) 

Using the value of tq determined in Sec. Ill CI we find that the variation in M^^. over the range of (3 values simulated 
on the LCP is up to 2% for the HISQ/tree action. We neglect the systematic effect introduced by this variation in 
the rest of the paper as it is of the same order as the statistical errors. The LCP for the asqtad action, however, 
corresponds to a strange quark mass that is about 20% heavier than the physical value. We will comment on how we 
account for this deviation from the physical value in Sec. |Vl 

Lattice estimates of hadron masses should agree with the corresponding experimental values in the continuum 
limit ^; however, at the finite lattice spacings used in thermodynamic calculations there are significant discretization 
errors. In staggered formulations, all physical states have taste partners with heavier masses that become degenerate 
only in the continuum limit. The breaking of the taste symmetry, therefore, introduces additional discretization 
errors, in particular, in thermodynamic observables at low temperatures where the degrees of freedom are hadrons. 
These artifacts have been observed in the deviations between lattice results and the hadron resonance gas model in 
the trace anomaly [l^ and in fluctuations of conserved charges [1^. In this subsection, we will quantify these taste 
symmetry violations in the asqtad, stout and HISQ/tree actions and show that they are the smallest in the HISQ 
action |46| . 

To discuss the effects of taste symmetry violations, we analyze all sixteen pseudoscalar mesons that result from 
this four-fold doubling, and are classified into eight multiplets with degenerate masses. These are labeled by their 
taste index F^ = 75, 7075, 7i75, 7oj 7ij lilo, lilji 1 [S^l. There is only one Goldstone boson, F^ = 75, that is 
massless in the chiral limit and the masses of the other fifteen pseudoscalar mesons vanish only in the chiral and 
continuum limits. The difference in the squared mass of the non-Goldstone and Goldstone states, — Afg, is the 
largest amongst mesons and their correlators have the best statistical signal; therefore, it is a good measure of taste 
symmetry violations. For different staggered actions, these violations, while formally of order a"a^, are large as 
discussed below ^ . 

The taste splittings, M^ — Mq, have been studied in detail for the asqtad and p4 actions [iO,!!^!!^!. The conclusion 
is that at a given /3 they are, to a good approximation, independent of the quark mass. Therefore, for the HISQ/tree 
action we calculate them on 16'^ x 32 lattices with mi/mg = 0.2 (see Table |V| . and on four 32'* ensembles with 
mi = 0.05ms for sea quarks and mi/mg = 0.2 for valence quarks. (The lattice parameters for these ensembles at 
/? = 6.664, 6.8, 6.95, and 7.15 arc given in Table IVll) The corresponding results, plotted in Fig. SI show the expected 
a^a^ scaling, similar to that observed previously with the HISQ action in the quenched approximation and in 
full QCD calculations with four flavors [i^H^. In this analysis, following Ref. [11], we use ay(q = 3.33/a) from 



^ For the nucleon and Q-baryon this has been demonstrated in Ref. [4l| . 

^ For the unimproved staggered fermion action as well as for the stout action the quadratic pseudoscalar meson splittings arc formally of 
order asO?, while for the asqtad and the HISQ actions they are of order a^a^. Projecting fat links to U{3) reduces the coefficients of 
the a"a^ taste violating terms. 
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Figure 4: The splitting M"^ — Mq of pseudoscalar meson multiplets calculated with the HISQ/tree and stout actions as a 
function of a%a^ (left). The right panel shows the RMS pion mass with Mg = 140 MeV as a function of the lattice spacing for 
the asqtad, stout and HISQ/tree actions. The band for the asqtad and stout actions shows the variation due to removing the 
fourth point at the largest a in the fit. These fits become unreliable for a>0.16 fm and are, therefore, truncated at a = 0.16 
fm. The vertical arrows indicate the lattice spacing corresponding to T ~ 160 MeV for A^^ = 6, 8 and 12. 



the potential as an estimate of as. Linear fits in ala? to the four points at the smallest lattice spacings shown in 
Fig. Ul^left) extrapolate to zero within errors in the eontinuum limit. The data also show the expected approximate 
degeneracies between the multiplets that are related by the interchange 7^ to 70 in the definition of as predicted 
by staggered chiral perturbation theory [62| . 

The splittings for the stout action, taken from Ref. [2^, for = 7^75 and 7^7^ are also shown in Fig. U] with open 
symbols. We find that they are larger than those with the HISQ/tree action for comparable lattice spacings. 

To further quantify the magnitude of taste-symmetry violations, we define, in MeV, the root mean square (RMS) 
pion mass as 

M^^'S = ^ ^ + + 3A/2 + 3M2 + 3A/2 + 3A/2 + + A/f) , (5) 

and plot the data in Fig. HTright) with Mq tuned to 140 MeV. The data for the asqtad and stout actions were taken 
from Ref. [s^ and Ref. [23j, respectively. As expected, the RMS pion mass is the largest for the asqtad action and 
smallest for the HISQ/tree action. However, for lattice spacing a ^ 0.104 fm, which corresponds to the transition 
region for Nt = 12, the RMS pion mass becomes comparable for the asqtad and stout actions. The deviations from 
the physical mass, Af^ = 140 MeV, become significant above a = 0.08 fm even for the HISQ/tree action. For the 
lattice spacings ^ 0.156 fm (a ^ 0.206 fm), corresponding to the transition region on — 8 {Nr = 6) lattices, the 
RMS mass is a factor of two (three) larger. 

Next, we analyze the HISQ/tree data for pion and kaon decay constants, given in Appendix [Cl for mi/nis = 0.05. We 
also analyze the fictitious rjss meson following Ref. [soj . In Fig.[51 we show our results in units of tq and ri determined 
in Sec. Ill CI as a function of the lattice spacing together with a continuum extrapolation assuming linear dependence 
on a^. We vary the range of the lattice spacings used in the fit and take the spread in the extrapolated values as an 
estimate of the systematic errors. These extrapolated values agree with the experimental results within our estimated 
errors (statistical and systematic errors are added in quadrature) as also shown in Fig. [S] This consistency justifies 
having used the continuum extrapolated value of /tt?'! from Ref. (58| to convert ri to physical units as discussed in 
Sec. Ill CI The deviation from the continuum value in the region of the lattice spacings corresponding to our finite 
temperature calculations is less than 8% for all the decay constants. We use these data to set the Jk scale and analyze 
thermodynamic quantities in terms of it and to make a direct comparison with the stout action data [23 - l23 | . 

Finally, in Fig. |6]we show the masses of cf) and K* mesons given in Appendix [Cl as a function of the lattice spacing. 
(The rho meson correlators are very noisy, so we do not present data for the rho mass.) Using extrapolations linear 
in we obtain continuum estimates, and by varying the fit interval, we estimate the systematic errors and add these 
to the statistical errors in quadrature. These estimates, in units of rg and ri, are plotted with the star symbol in 
Fig. |6l The experimental values along with error estimates are shown as horizontal bands and agree with lattice 
estimates, thereby providing an independent check of the scale setting procedure. The slope of these fits indicates 
that discretization errors are small and confirms the findings in [46| that taste symmetry violations are much smaller 
in the HISQ/tree action compared to those in the asqtad action. For the range of lattice spacings relevant for the 
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Figure 5; The decay constants of rjss, K and n mesons with the HISQ/tree action at mi — 0.05ms measured in units of ro 
(left) and ri (right) are shown as a function of the lattice spacings. The black points along the y-axis are the result of a linear 
extrapolation to the continuum limit. The experimental results are shown as horizontal bands along the y-axis with the width 
corresponding to the error in the determination of ro and ri, respectively. We use the HPQCD estimate for fr,^g [1^ as the 
continuum value. 
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Figure 6: The masses of the (fi and K* mesons with the HISQ/tree action at mi — 0.05ms measured in units of ro (left) 
and ri (right) are shown as a function of the lattice spacing. The lines show a linear continuum extrapolation and the black 
crosses denote the extrapolated values. The experimental results are shown as horizontal bands along the j/-axis with the width 
corresponding to the error in the determination of ro and ri, respectively. 



finite temperature transition region on Nr = 6-12 lattices, the discretization errors in the vector meson masses are 
less than 5%. 



III. UNIVERSAL SCALING IN THE CHIRAL LIMIT AND THE QCD PHASE TRANSITION 

In the limit of vanishing light quark masses and for sufficiently large values of the strange quark mass, QCD is 
expected to undergo a second order phase transition belonging to the universality class of three dimensional 0(4) 
symmetric spin models 0] . Although there remains the possibility that a fluctuation-induced first order transition may 
appear at (very) small values of the quark mass, it seems that the QCD transition for physical values of the strange 
quark mass is, indeed, second order when the light quark masses arc reduced to zero. An additional complication in 
the analysis of the chiral phase transition in lattice calculations arises from the fact that the exact 0(4) symmetry 
is difficult to implement at nonzero values of the lattice spacing. Staggered fermions realize only a remnant of this 
symmetry; the staggered fermion action has a global 0(2) symmetry. The restoration of this symmetry at high 
temperatures is signaled by rapid changes in thermodynamic observables or peaks in response functions, which define 
pseudocritical temperatures. For these observables to be reliable indicators for the QCD transition, which becomes a 
true phase transition only in the chiral limit, one must select observables which, in the chiral limit, are dominated by 
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contributions arising from the singular part of the QCD partition function Z{V,T), or more precisely from the free 
energy density, / = — TV~^ In Z{V,T). A recent analysis of scaling properties of the chiral condensate, performed 
with the p4 action on coarse lattices, showed that critical behavior in the vicinity of the chiral phase transition is well 
described by 0{N) scaling relations [6^ which give a good description even in the physical quark mass regime. 

In the vicinity of the chiral phase transition, the free energy density may be expressed as a sum of a singular and 
a regular part, 

T 

f = -—\nZ = fsing{t,h) + freg{T,nii,ms) . (6) 

Here t and h are dimcnsionless couplings that control deviations from criticality. They are related to the temperature 
T and the light quark mass mi , which couples to the symmetry breaking (magnetic) field, as 

to "0 rus 

where denotes the chiral phase transition temperature, i.e., the transition temperature at _ff = 0. The scaling 
variables t, h are normalized by two parameters to and ho, which arc unique to QCD and similar to the low energy 
constants in the chiral Lagrangian. These need to be determined together with T^. In the continuum limit, all three 
parameters are uniquely defined, but depend on the value of the strange quark mass. 

The singular contribution to the free energy density is a homogeneous function of the two variables t and h. Its 
invariance under scale transformations can be used to express it in terms of a single scaling variable 

where /3 and S are the critical exponents of the 0{N) universality class and zq — to/h]J^^ . Thus, the dimcnsionless 
free energy density / = f /T* can be written as 

fiT,mi,ms) = h'+'/\fsiz) + friT,H,ms) , (9) 

where the regular term gives rise to scaling violations. This regular term can be expanded in a Taylor series around 
{t,h) = (0,0). In all subsequent discussions, we analyze the data keeping iris in Eq. ^ fixed at the physical value 
along the LCP. Therefore, the dependence on nis will, henceforth, be dropped. 

We also note that the reduced temperature t may depend on other couplings in the QCD Lagrangian which do not 
explicitly break chiral symmetry. In particular, it depends on light and strange quark chemical potentials fiq, which 
in leading order enter only quadratically, 

/ rp rjiQ , ,, ^ 9. ... ,. \ 

(10) 

Derivatives of the partition function with respect to are used to define the quark number susceptibilities. 

The above scaling form of the free energy density is the starting point of a discussion of scaling properties of most 
observables used to characterize the QCD phase transition. We will use this scaling Ansatz to test to what extent 
various thermodynamic quantities remain sensitive to universal features of the chiral phase transition at nonzero 
quark masses when chiral symmetry is explicitly broken and the singular behavior is replaced by a rapid crossover 
characterized by pseudocritical temperatures (which we label T^) rather than a critical temperature. 

A good probe of the chiral behavior is the 2-flavor light quark chiral condensate 

Following the notation of Ref. [G^I , we introduce the dimcnsionless order parameter Mi 




b, 

Mb . ^^^^^^^ . (12) 



Multiplication by the strange quark mass removes the need for multiplicative renormalization constants; however, Mb 
does require additive renormalization. For a scaling analysis in h at a fixed value of the cutoff, this constant plays no 
role. Near T°, Aff, is given by a scaling function /g(-z) 

Mt{T, H) - h^/'fait/h^'f"') + lM,reg{T, H) , (13) 
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and a regular function fM,reg{T,H) that gives rise to scaling violations. We consider only the leading order Taylor 
expansion of fM.regiT, H) in H and quadratic in t, 



fALregiT,H) = at{T)H 

= («o + ai^^+a2(^^^j jH (14) 

with parameters cq, ai and 02 to be determined. The singular function fa is well studied in three dimensional spin 
models and has been parametrized for the 0(2) and 0(4) symmetry groups [GSl - lssj . Also, the exponents /3, 7, 5 and 
v used here are taken from Table 2 in Ref. |68j . 

Response functions, derived from the light quark chiral condensate, are sensitive to critical behavior in the chiral 
limit. In particular, the derivative of with respect to the quark masses gives the chiral susceptibility 

xm,i = ^^i^^yi'^' ■ (15) 

The scaling behavior of the light quark susceptibility, using Eq. is 



f ^ , 1/5-1 f , dfM,reg{T,H) 



Xm,l 

with /^(z)=:i[/G(z)-|/^(z)]. (16) 

The function has a maximum at some value of the scaling variable z = Zp. For small values of h this defines the 
location of the pseudocritical temperature as the maximum in the scaling function fciz)- Approaching the critical 
point along h with z fixed, e.g., z = Q 01 z = Zp, Xm,i diverges in the chiral limit as 



Similarly, the mixed susceptibility 



Xm,i^m/ . (17) 



also has a peak at some pseudocritical temperature and diverges in the chiral limit as 

(/3-l)/,35 

Xt,i~S '"^ . (19) 

One can calculate xt,i either by taking the derivative of {■ipip) with respect to T or by taking the second derivative 
with respect to /ij, i.e., by calculating the coefficient of the second order Taylor expansion for the chiral condensate 
as a function of fj,i /T [63 • The derivative of ('0V') with respect to T is the expectation value of the chiral condensate 
times the energy density, which is difficult to calculate in lattice simulations, as additional information on temperature 
derivatives of temporal and spatial cutoff parameters is needed. Taylor expansion coefficients, on the other hand, are 
well defined and have been calculated previously, although their calculation is computationally intensive. This mixed 
susceptibility has been used to determine the curvature of the chiral transition line for small values of the baryon 
chemical potential [69| . 

Other thermodynamic observables analyzed in this paper are the light and strange quark number susceptibilities 
defined as 

Xg 1 dHnZ 

T2 - vT^^{^iJTY ' ^ • ^^"^ 

These are also sensitive to the singular part of the free energy since the reduced temperature t depends on the quark 
chemical potentials as indicated in Eq. (jlOp . However, unlike the temperature derivative of the chiral condensate, i.e., 
the mixed susceptibility Xt,l^ the temperature derivative of the light quark number susceptibility does not diverge in 
the chiral limit. Its slope at T° is given by 
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and has the contribution Cr from the regular part of the free energy, while its variation with temperature is controlled 
by the singular part. The critical exponent a is negative for QCD since the chiral transition is expected to belong to 
the universality class of three-dimensional 0{N) models. In short, while Xq is sensitive to the critical behavior, it does 
not diverge in the thermodynamic limit. Consequently, it has been extremely difficult to extract reliable information 
on T° or Tc from scaling fits to Xq- Even in high statistics 0{N) model calculations [t^ the structure of the subleading 
term in Eq. pip could only be determined after using results for the dominant contribution Cr extracted from other 
observables. We, therefore, consider quark number susceptibilities as a good indicator of the transition in QCD, but 
not useful for extracting precise values for the associated pseudocritical temperature. 
Finally, we consider the expectation value of the Polyakov loop L, 

1 

L{x) = -Tyl[Uo{xo,x) , (22) 

Xo = l 

which is the large distance limit of the static quark correlation function, 

= Jim (i(O)L^(f)) . (23) 

L = {L{x)) is a good order parameter for deconfinement in the limit of infinitely heavy quarks. In that limit, it can 
be related to the singular structure of the partition function of the pure gauge theory and can be introduced as a 
symmetry breaking field in the action. In QCD with light quarks, in particular in the chiral limit, L is no longer an 
order parameter due to the explicit breaking of the Z(Z) center symmetry by the quark action. It does not vanish for 
T <T^, but is determined by the value of the free energy of a static quark Fq in the confined phase. This free energy 
can be well approximated by the binding energy of the lightest static-light meson, which is of order Aqcd ■ Similarly, 
Tc ^ Aqcd', consequently L ~ exp(— i^g/T) ~ 1/e is not small in the confined phase. The data for QCD with 
light quarks show that L varies significantly with temperature in the transition region, reflecting the rapid change in 
screening properties of an external color charge. Thus, the Polyakov loop is sensitive to the transition but it has no 
demonstrated relation to the singular part of the QCD partition function. We, therefore, do not use it to determine 
an associated pseudocritical temperature. 



IV. CHIRAL OBSERVABLES 



In this section, we present results for observables related to chiral symmetry restoration at finite temperatures and 
discuss the cutoff dependence of these quantities. To set the normalization of different quantities we express them in 
terms of the staggered fermion matrix Dq = niq ■ 1 + D with q = l,s, as in Ref. [l^j. In what follows, {'4'ij)q,T will 
denote the one-flavor chiral condensate, i.e., 

W),,. = J^^Tr(i?-i), q = l,s, (24) 

where the subscript x = t and a; = will denote the expectation value at finite and zero temperature, respectively. 
The chiral susceptibility defined in Sec. Illll is the sum of connected and disconnected Feynman diagrams defined as 

Xm,l(J- ) = 2 = XI, disc + XI, con , (25) 

onii 

Xra^sx-^ I — Xs,disc ~r As, con ; 



drria 



with 

Xq 



-[{{TrD-^)^) - {TTD-^f) , (27) 



and 



Xq,con ^ -'^TtY,{D;\x,0)D-\0,x)) , q = l,s. (28) 



Here ti/ = 2 for light quark susceptibilities, and = 1 for the strange quark susceptibilities. The disconnected part 
of the light quark susceptibility describes the fluctuations in the light quark condensate and is directly analogous to 
the fluctuations in the order parameter of an 0{N) spin model. The second term (xq.con) arises from the explicit 
quark mass dependence of the chiral condensate and is the expectation value of the volume integral of the correlation 
function of the (isovector) scalar operator ipip. 
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Figure 7: The subtracted chiral condensate for the asqtad and HISQ/tree actions with mi = ms/20 is compared with the 
continuum extrapolated stout action results [S^] (left panel). The temperature T is converted into physical units using ri in 
the left panel and /k in the right. We find that the data collapse into a narrow band when fx is used to set the scale. The 
black diamonds in the right panel show HISQ/tree results for Nt = 8 lattices after an interpolation to the physical light quark 
mass using the mi/ms = 0.05 and 0.025 data. 



A. The chiral condensate 



The chiral condensate (■0V') requires both multiplicative and additive renormalizations at finite quark masses. The 
leading additive renormalization is proportional to (m^/a^).'* To remove these UV divergences, we consider the 
subtracted chiral condensate introduced in Ref. [l9| . 

Our results for the HISQ/tree and asqtad actions at mi = 0.05ms are shown in Fig. [Tl^left) and compared to the 
continuum estimate obtained with the stout action [2^ . The temperature scale is set using tq and ri as discussed 
in Section fll CI and Appendix iBl The asqtad results obtained on Nr — 8 lattices deviate significantly from the stout 
results as observed previously [201. The new data show that these differences are much smaller for A^,- = 12 ensembles. 
More important, the discretization effects and the differences from the stout continuum results arc much smaller for 
the HISQ/tree data. 

In Fig. [fright), we analyze the data for A/^^ using the kaon decay constant Jk to set the lattice scale. For 
the HISQ/tree action, we use the values of Jk discussed in Sec. IlICi while for the asqtad action, we use fx from 
staggered chiral fits (see the discussion in appendix B). We note that for mi/ms = 1/20 all the data obtained with 
the HISQ/tree and asqtad actions on different Nr lattices collapse into one curve, indicating that A/^^ and /k have 
similar discretization errors. The remaining difference between the stout and our estimates, as shown next, is due 
to the difference in the quark masses — calculations with tire stout action were done with mi = 0.037ms whereas our 
calculations correspond to m; = O.OSm^. 

For a direct comparison with stout results, we extrapolate our HISQ/tree data in the light quark mass. This 
requires estimating the quark mass dependence of the chiral condensate at both zero and non-zero temperatures. For 
the r = data, we perform a linear extrapolation in the quark mass using the HISQ/tree lattices at m; = O.OSm^ 
and 0.20m.,. For the non-zero T data, we use the 0{N) scaling analysis, described in Sec. |Vl which gives a good 
description of the quark mass dependence in the temperature interval 150MeV < T < 200 McV. The resulting Nr ~ 8 
HISQ/tree estimates at the physical quark mass are shown in Fig. [TKright) as black diamonds and agree with the 
stout action results [l^-HJ plotted using green triangles. 

We can also remove the multiplicative renormalization factor in the chiral condensate by considering the renormal- 
ization group invariant quantity rfms{'tptlj)i, where is the strange quark mass. The additive divergences can be 



* There is also a logarithmie divergenee proportional to m^, which we neglect. 
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Figure 8: The renormalized chiral condensate for the HISQ/tree action with mi/nis = 0.05 is compared to the stout data. 
In the right panel, we show the renormalized strange quark condensate for the HISQ/tree action. The temperature scale 
in both figures is set using r\. The black diamonds in the left panel show the Nt ~ 8 HISQ/tree estimates using the /k scale 
and after an interpolation to the physical light quark mass mi/nis — 0.037 as discussed in the text. 



removed by subtracting the zero temperature analogue, i.e., we consider the quantity 

Af = d + 2m,rt{{{^i^)i,r - W)/,o) • (30) 

Note that Af is very similar to the renormalized chiral condensate introduced in Ref. j2^ . but differs by 

the factor {nii/ms)/{rfm^) and d. A natural choice for d is the value of the chiral condensate in the zero light 
quark mass limit times nigrf. In this limit Af should vanish above the critical temperature. To estimate d, we use 
the zero temperature estimate {r{jip)i{MS, fi = 2GeV) ~ 242 (9) (j]^?) (4) MeV^ determined in the chiral limit using 
SU{2) staggered chiral perturbation theory by the MILC collaboration [4l| and the corresponding strange quark mass 
m^ifl = 2GeV) ^ 88(5) MeV. We get d = 0.0232244. 

We show Af for the HISQ/tree action and the stout continuum results in Fig. [SJleft)^. To compare with the stout 
continuum results, we need to extrapolate the HISQ/tree data both to the continuum limit and to the physical quark 
mass. To perform the continuum extrapolation we convert Af to the Jk scale in which discretization errors, as 
already noted for A/^g, are small. We then interpolate these Nr ~ 8 data at mi/nis = 0.05 and 0.025 to the physical 
quark mass mi/ms = 0.037. These estimates of the continuum HISQ/tree Af are shown in Fig. |8l[left) as black 
diamonds and arc in agreement with the stout results (green triangles) |24| . 

Lastly, in Fig. [S^right), we show the subtracted renormalization group invariant quantity. Af , which is related 
to the chiral symmetry restoration in the strange quark sector. We find a significant difference in the temperature 
dependence between Af and Af , with the latter showing a gradual decrease rather than a crossover behavior. 



B. The chiral susceptibility 



As discussed in Sec. lIIIl the chiral susceptibility Xm.i is a good probe of the chiral transition in QCD as it is sensitive 
to the singular part of the free energy density. It diverges in the chiral limit, and the location of its maximum at 
nonzero values of the quark mass defines a pseudocritical temperature Tc that approaches the chiral phase transition 
temperature T° as m/ — 0. 

For sufficiently small quark masses, the chiral susceptibility is dominated by the disconnected part, therefore, Tc can 
also be defined as the location of the peak in the disconnected chiral susceptibility defined in Eq. (P7|) . As we will show 
later, Xq,disc does not exhibit an additive ultraviolet divergence but does require a multiplicative renormalization ^ . 



^ We multiply the stout results by {nis/mi) = 27.3 and by rfm'^ = 0.0022275. For the latter factor, we use the physical pion mass and 

the value of ri determined in [5ql and discussed in Sec. Ill CI 
® It is easy to see that at leading order in perturbation theory, i.e., in the free theory, the disconnected chiral susceptibility vanishes and 

thus is non-divergent. Our numerical results at zero temperature do not indicate any quadratic divergences in the disconnected chiral 

susceptibility, but logarithmic divergences are possible. 
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Figure 9: Value of the bare light quark mass, in MeV using ri to set the scale, on the line of constant physics for the asqtad 
action vs. the lattice gauge coupling /?. The right hand part of the figure shows the change of the renormalization constants 
with P, I.e., with the cutoff, relative to the arbitrarily chosen renormalization point /3 — 6.65. 



1. Disconnected chiral susceptibility 



The multiplicative renormalization factors for the chiral condensate and the chiral susceptibility can be deduced 
from an analysis of the line of constant physics for the light quark masses, m;(/3). The values of the quark mass for 
the asqtad action, converted to physical units using ri, arc shown in Fig. [9l^left). The variation with /3 gives the scale 
dependent renormalization of the quark mass (its reciprocal is the renormalization factor for the chiral condensate). 
What mi{p) does not fix is the renormalization scale, which we choose to be ro/a = 3.5 (equivalently ri/a = 2.37 
or a = 0.134 fm), and the "scheme", which we choose to be the asqtad action. For the asqtad action, this scale 
corresponds to the coupling /3 = 6.65 which is halfway between the peaks in the chiral susceptibility on N^- = 8 
and 12 lattices. This specification, Zm{asqtad) = 1 at ro/a = 3.5, is equivalent to choosing, for a given action, the 
renormalization scale A which controls the variation of Zm with coupling /3 as shown in Fig. [9][riglit) for the asqtad 
action. 

A similar calculation of Zm is performed for the p4 and HISQ/tree actions. It is important to note that choosing 
the same reference point ro/a = 3.5 and calculating Zm(/3) for each of the actions leaves undetermined a relative 
renormalization factor between the actions, i.e., the relation between the corresponding A's of the different schemes. 
This relative factor between any two actions is also calculable and given by the ratio of the (bare) quark mass along 
the physical LCP at ro/a = 3.5. At this scale our data give 



m(asqtad) 
m (HISQ/tree) 



0.97828. (31) 



Recall, however, that along the LCP the quark masses, and therefore mi, for the asqtad action are about 20% 
heavier than the physical values. Noting that the lattice scale at a given (3 is set using a quark-mass independent 
procedure, we correct TO(asqtad) by the factor (A/^ro)^|HisQ/troc/(A'j'^7r''o)^|asqtad- Then, at ro/a = 3.5 

m(asqtad) Z„, (asqtad) , , 

^ ^ ' 0.782 i.e. 7uTQo/ ■ = 1.2786 . (32) 



m(HISQ/tree) ' ' Z„ (HISQ/tree) 

Given Zm{l3) we get Z^^, = Zs = 1/^m and Z^ = l/Z^,- A similar calculation of Z,„ has been carried out for the p4 
action. 

The systcmatics of the quark mass and cutoff dependence of the disconnected part of the chiral susceptibility is 
analyzed in more detail for the p4 and asqtad actions in Fig. 1101 The data show a rapid rise in xi,disc/T'^ with 
decreasing quark mass at low temperatures and in the transition region. This mass dependence can be traced back to 
the leading thermal correction to the chiral condensate. At finite temperature and for sufficiently small quark masses, 
the chiral order parameter can be understood in terms of the 3-dimensional 0{N) models. A dimensional reduction is 
applicable because the Goldstone modes are light in this region. Based on the 0{N) model analysis, the quark mass 
dependence of the chiral condensate is expected to have the form [7l|-[74| 



(V^V)i(r,m,) = (i^^>i(0) + C2(T)V^+... , 



(33) 
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Figure 10: The disconnected part of the chiral susceptibility, including multiplicative renormalization, calculated on Nt — 8 
lattices for the p4 (left) and asqtad (right) actions at three light quark masses. The figure on the right also shows asqtad data 
from Nt = 6 lattices as open symbols. 

as has been confirmed in numerical simulations with the p4 action on Nr ~ 4 lattices (63 |. Consequently for T < T^, 

— 1/2 

there is a singularity in the chiral susceptibility in the limit of zero quark mass which explains the rise in 

\l^disc/T . 

A second feature of the data is shown in Fig. [TOl right) which compares data for the asqtad action on lattices of 
different at mi/ms = 0.2 and 0.1. Open (filled) symbols denote data on Nr = 6 {Nr = 8) lattices. The data show 
a shift towards smaller temperature values of both the peak and the rapidly dropping high temperature part when the 
lattice spacing is reduced. The data also show that the variation of the shape of the susceptibility above the peak is 
weakly dependent on the quark mass. This is expected as xi is the derivative of the chiral condensate with respect to 
the mass which, as shown in Fig. [8l is almost linear in the quark mass in this temperature regime. Thirdly, the data 
in Fig. [TUTright) show that the height of the multiplicatively renormalized disconnected chiral susceptibility at fixed 
mi/iris is similar for Nr = 6 and Nr = 8 lattices. This lack of increase in height with Nr supports the hypothesis 
that there are no remaining additive divergent contributions in the disconnected part of the chiral susceptibility. 

In Fig. Ill[ we compare, for mi/nis = 0.05, the disconnected part of the chiral susceptibility including the multi- 
plicative renormalization factor Z-,^. We note three features in the data. First, the variation in the position of the peak 
for the asqtad action is larger between Nr = 8 and 12 than for the HISQ/tree action between Nr ~ 6 and 8. Second, 
the peak height increases for the HISQ/tree data and decreases for the asqtad data with Nr- Lastly, the agreement 
in the location of the peak for the two actions and the data above the peak is much better when /k is used to set 
the scale as shown in Fig. [TlTright). Note that the peak height for the two actions is not expected to match since the 
quark masses on the LCP for the asqtad data are about 20% heavier than for the HISQ/tree data. 

2. Connected chiral susceptibilities 

The connected part of the chiral susceptibility, Eq. (|28|) . is the volume integral of the scalar, flavor nonsinglet 
meson correlation function. At large distances, where its behavior is controlled by the lightest scalar, nonsinglet 
screening mass, the correlation function drops exponentially as this state has a mass gap, i.e., xi.con can diverge in 
the thermodynamic limit only if the finite temperature screening mass in this channel vanishes. This, in turn, would 
require the restoration of the Ua{^) symmetry, which is not expected at the QCD transition temperature. In fact, 
the scalar screening masses are known to develop a minimum at temperatures above, but close to, the transition 
temperature (TSj . One therefore expects that even in the chiral limit, the connected part of the chiral susceptibility 
will only exhibit a maximum above the chiral transition temperature. 

There are two subtle features of the connected part of the susceptibility calculated at nonzero lattice spacings that 
require further discussion. First, taste symmetry violations in staggered fermions introduce an additional divergence 
of the form / y/ml for T < T°. It also arises due to the long distance fluctuations of Goldstone pions, as explained 
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Figure 11: The disconnected part of the chiral susceptibility for the asqtad and HISQ/tree actions, including the multiplicative 
renormalization constant discussed in the text, is shown for mi = ma/20 and different A'^^- In the right panel, the same data 
are plotted using /x to set the scale. Plotted this way they show much smaller variation with A^'^. 



in Eq. (33), however, unlike the divergence in the disconnected part which is physical, this term is proportional to 
the O(a^) taste breaking. Note that in the two- flavor theory there are no such divergences due to Goldstone modes 
in the continuum limit jT^, [z3|- Thus, wc expect to observe a strong quark mass dependence at low temperatures in 
Xi,con- Second, there is a large reduction in the Ua{^) symmetry breaking in the transition region, consequently there 
will be a significant quark mass dependence of scalar screening masses and of xi,con- 

We have calculated xi.con for the p4, asqtad and HISQ/tree actions. Results at different light quark masses from 
Nt = 8 lattices are shown in Fig. [T^] for the p4 and asqtad actions with the multiplicative renormalization performed 
in the same way as for xi,disc- A strong dependence on the quark mass is seen in both the p4 and asqtad data. 
This, as conjectured above, is due to a combination of the artifacts that are due to taste symmetry breaking and the 
variations with temperature of the scalar, flavor nonsinglet screening mass at and above the crossover temperatures. 

In Fig. [lH we show the connected chiral susceptibility for the asqtad and HISQ/tree actions at fixed mi ~ 0.05ms 
for different Nr- In Sec. IIV A} we noted the presence of an additive quadratic divergence, proportional to mq/a^, 
in the chiral condensate, which will give rise to a mass-independent quadratic divergence in the chiral susceptibility. 
We find that the absolute value of the data grows with Nt as expected. Since this divergent contribution is the 
same for light and strange susceptibilities, it can be eliminated by constructing the difference xi,con — '2Xs,con- The 
resulting data arc shown in Fig. Upright), and wc find that the peak occurs at slightly higher T as compared to the 
disconnected chiral susceptibility shown in Fig. [TT] Also, wc find that the height of the peak decreases with Nr and 
the position of the peak is shifted to smaller temperatures on decreasing the lattice spacing, which is most evident 
when comparing the Nr = 8 and Nr — 12 asqtad data. 



3. Renormalized two-flavor chiral susceptibility 

Lastly, we compare our estimates for the two-flavor chiral susceptibility, deflned in Eqs. (P7)) and with results 
obtained with the stout action To remove the additive ultraviolet divergence discussed above, we now subtract 
the zero temperature rather than the strange quark chiral susceptibility. Furthermore, to get rid of the multiplicative 
renormalization, this combination is multiplied by rnj., i.e., the following quantity is considered 

= ^ ix^AT) - X^AT = 0)) . (34) 

This construct has the advantage of being renormalization group invariant and, unlike the definition in Ref. (22j . it 
does not vanish in the chiral limit. In Fig. IWleft). we show data for the stout, HISQ/tree and asqtad actions with 
the scale set by n. The stout data have been taken from Ref. [22| and multiplied by {•nis/mif' = (27.3)^ to conform 
to Eq. [23.[23j. The HISQ/tree results on Nr = 12 lattices are not shown as the corresponding zero temperature 
calculations are not yet complete. We find that the large difference between the continuum stout and Nr = 8 asqtad 
results is significantly reduced by Nr = 12. Second, the cutoff dependence for the HISQ/tree data is much smaller 
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Figure 12: The connected part of the chiral susceptibihty for the p4 (left) and asqtad (right) actions for different quark masses 
on Nr = 8 lattices. 
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Figure 13: The connected part of the chiral susceptibility for the asqtad and HISQ/tree actions on the LCP defined by 
mi — 0.05ms. In the right hand figure, we show the difference of the light and strange quark connected susceptibilities in which 
the divergent additive artifact cancels. 



than for the asqtad data. A similar behavior was also observed in the case of the chiral condensate, as discussed in 
Sec. WJ\ 

The cutoff dependence between the HISQ/tree, asqtad and continuum stout data is significantly reduced when fx 
is used to set the scale as shown in Fig. [Tif right). The change in the HISQ/tree data is small as the scales determined 
from ri and /k are similar. The difference in the scales from the two observables is larger for the asqtad action at 
these lattice spacings and using /k mostly shifts the Nr ~ 8 data. The difference in the position of the peak between 
the three actions also decreases, whereas the height of the peak and the value in the low temperature region show 
significant differences between the stout and HISQ/tree (or asqtad) actions. Since Xfl(T)/T^ is a renormalization 
group invariant quantity, the only reason for the difference should be the different values of the light quark mass: 
mi/ms ~ 0.05 for the HISQ/tree estimates versus 0.037 for the stout data. 0{N) scaling, discussed in Sec. IIII[ 
suggests that the peak height should scale as h^/^~^ ^ m^'^'^. Applying this factor to the stout data reduces the 
peak height from ~ 40 to '-^ 31.5. Similarly, two corrections need to be applied to the asqtad data. First, a factor 
of 1/1.44 to undo the multiplication by a heavier and the second, a multiplication by a factor of 1.2 to scale the 
susceptibility to the common light quark mass. After making these adjustments to normalize all three data sets to 
mi/ms = 0.05, we find that the HISQ/tree, asqtad and stout data are consistent. 
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Figure 14: The renormalized two-flavor chiral susceptibility xr for the asqtad and HISQ/tree actions obtained at mi — 0.05ms 
and compared with the stout action results [2^. The temperature scale is set using ri {fx) in the left (right) panels. 

V. 0(JV) SCALING AND THE CHIRAL TRANSITION TEMPERATURE 

A. The transition temperature using the p4 action 

In this section, wc use the universal properties of the chiral transition to define the transition temperature and 
its quark mass dependence for sufficiently small quark masses, as discussed in Sec. IIIIl The scaling analysis of the 
chiral condensate leads to a parameter free prediction for the shape and magnitude of the chiral susceptibility. In the 
vicinity of the chiral limit, the peak in the chiral susceptibility corresponds to the peak in the scaling function /^(z) 
and the quark mass dependence of the pseudocritical temperature Tc is controlled entirely by the universal 0{N) 
scaling behavior. Keeping just the leading term proportional to ai in the regular part, the position of the peak in 
Xm,! is determined from Eq. p6|) using 



dT V T4 



(35) 



which, for zero scaling violation term, i.e., ai = 0, gives the position of the peak in the scaling function at z ~ Zp 
(see Sec. IIIip . The strange quark mass on the left hand side is included only for consistency as the derivative is taken 
keeping it constant. For small light quark masses, we can expand fx{z) around Zp: 



In this approximation, the location of the maximum in the chiral susceptibility varies as 

^~ " 2Ap 

and the variation of the pseudocritical temperature as a function of the quark mass is given by 



(36) 



(37) 



ai 



2ApZpZohQ 



i/s 



I 2ApZpZo^-P 



(38) 



Recall that is the transition temperature in the chiral limit. Thus, to determine the pseudocritical temperatures 
Tc{H), we need to perform fits to the chiral condensate Mb, defined in Eqs. and to determine the parameters 
T°, zq, to, Qq, Oi and 02 in the scaling and regular terms. Theoretically, one expects the 0(4) Ansatz to describe the 
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Figure 15: Estimates of pseudocritical temperature determined from the peak in the disconnected susceptibility with the p4 
action on lattices of temporal extent Nt = 4 and 8 (left). Results for rris/mi < 0.05 will be presented elsewhere. Curves 
show results obtained from the 0(2) scaling fits to the chiral condensate without (blue) and with (red) scaling violation terms 
included. The right hand figure shows the corresponding analysis for the asqtad action. 



critical behavior in the continuum limit; however, 0{2) scaling is more appropriate for calculations with staggered 
fermions at non-zero values of the cutoff. Fits are, therefore, performed using both the 0{2) and 0(4) scaling functions 
in all cases. 

Performing universal 0{N) scaling fits to extract reliable estimates of the chiral transition temperature as a function 
of the light quark masses requires making choices for the range of data points (temperature values) to include and the 
number of terms needed to model the regular part. To test these issues numerically, we use the extensive data with the 
p4 action on lattices with Nr = 4 and 8 and light quark masses down to mi/ms = 1/80 [13, [6^ [t^. In Fig. irsT left). 
we compare scaling results for Tc{H) as a function of the light quark mass, H = mi/ms, with direct determination 
of the peak in the disconnected part of the chiral susceptibility obtained using Gaussian and cubic fits to the data 
in the peak region. We find that the scaling fits reproduce the quark mass dependence of Tc{H) for m//TOs<l/20, 
which covers the physical point mi/rrig ~ 1/27. For larger quark masses, scaling violation terms have to be taken into 
account. The range of temperature values included in these fits (defining the scaling window) is 15-20 MeV. 

Having demonstrated that the scaling analysis works for the p4 action, we apply this approach to the asqtad and 
HISQ/trec actions even though the range of quark masses explored is smaller and the coverage of the transition region 
(values of T simulated) is not as dense, so a larger T range will be used. 



B. Scaling analysis for the asqtad and HISQ/tree action 

The scaling analysis of the HISQ/tree {Nt = 6, 8 and 12 lattices) and asqtad actions [Nr — 8 and 12 lattices) 
was performed in steps due to the limited number of mi/ms and temperature values simulated compared to the p4 
action discussed in Sec. IV Al In the case of the HISQ/tree action, having data at mi/nis = 1/20 and 1/40 on Nr ~ 6 
and Nr = 8 lattices allowed us to test the range of validity of the scaling fits in the region bracketing the physical 
light quark mass. The scaling Ansatz included both the singular and regular parts defined in Eqs. and ([55)1 in 
terms of the six parameters T^, zo, ^Oi QOi oi and 02. We find that a best fit to the rai/ms = 1/40 data for Mi, also 
fits the mi/nis = 1/20 data; in addition, the chiral susceptibility derived from this fit matches the measurements at 
mi/ms = 1/20. We, therefore, conclude that the scaling extends to mi/nis = 1/20. 

To get the value of Tc, we then carried out a simultaneous fit to the chiral condensate at the two values of mi/uis 
and its derivative, the total chiral susceptibility, at mi/ms = 1/20. (The calculation of the connected part of the 
susceptibility at mi/ms = 1/40 is not complete, so we could not include these data.) The range of temperature values 
selected for the three data sets was independently adjusted to minimize the x^/dof and, at the same time, include as 
many points as possible in the region of the peak in the susceptibility. We also varied the relative weight assigned to 
points along the two LCPs and for our final fits decreased the weight given to the heavier mi/rris = 1/20 points by a 
factor of sixteen in the calculation of the ■ (This was done because the scaling Ansatz with a truncated regular piece 
is expected to get progressively worse as the mass increases, whereas the statistical errors in Mb are roughly a factor 
of two smaller at the heavier mass.) In all these fits we find that while the height of the peak in the susceptibility 
is sensitive to what points are included in the fits and the relative weighting of the points along the two LCP, the 
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Figure 16: Scaling fits and data for the chiral condensate Mb calculated with the HISQ/tree action on lattices with temporal 
extent A',- = 6 (left) and the chiral susceptibility Xm,! (right). The data for Mh at mi/ms = 0.025 and for Mb and Xm,i at 
mi I ma = 0.05 are fit simultaneously using the 0(4) scaling Ansatz. The fits using the 0(2) Ansatz are similar. The points 
used in the scaling fits are plotted using open symbols. The dotted lines give the data scaled to the physical quark masses. 

location of the peak does not vary by more than 0.3 MeV for any reasonable set of choices. Repeating the analysis 
vifith just the mi /iris = 1/40 or the mi /iris = 1/20 data did not change the estimates of Tc significantly. We therefore 
conclude that mi /iris = 1/20 lies within the range of validity of our scaling Ansatz. Confirming that mi/mg = 1/20 
lies within the scaling window is important for our analysis because simulations on the Nt = 12 lattices have been 
done only at this LCP for both the HISQ/tree and asqtad actions. 

The most challenging part of determining the best fit was including a sufficient number of points below the location 
of the peak in Xm,i- This is evident from Figs. [TMTQI which show that the fits rapidly deviate from the Xm,i data 
below the location of the peak. In most cases, only two points below the peak could be included in the fits. 

The x^/dof of these fits are not good. For the 0(4) fits they are 96/12 and 30/12 for the asqtad data on A^,- = 8 
and 12 lattices; and 430/28, 125/26 and 70/14 for the iVr = 6, 8 and 12 HISQ/tree data, respectively. (The x^/dof 
for the 0(2) fits is about 20% larger for the A^^. = 6 and 12 HISQ/tree data and comparable for the others.) These 
large x^/dof reflect the fact that the statistical errors in Mf, arc small and do not include the systematic errors due 
to fluctuations in points along the LCP coming from less than perfect tuning of rrig (or T), and partly because we 
have used a truncated form for the regular part of the free energy. We attempted various truncations of the regular 
part and found that the variation in the location of the peak is negligible and insensitive to the number of parameters 
included in the fits. Note that the location of the peak and its scaling with mi/mg are the only quantities we extract 
from the fits. Overall, we find that the fits using the 0(4) scaling Ansatz are more stable with respect to variation of 
the fit range. We, therefore, use estimates from the 0(4) fits to obtain our final results for Tc which are about 2 MeV 
lower than those with the 0(2) Ansatz. Fits to the HISQ/tree data are shown in Figs. [121 [T71 and [TBI for Nr = 6, 8 
and 12, respectively. 

The above analysis using simultaneous fits to data for the order parameter and the total chiral susceptibility x 
was repeated for the asqtad action on Nr = 8 lattices with the two LCP at mi/ms = 1/20 and 1/10. We failed to 
find a fit that reproduced the data for the chiral susceptibility at the heavier mass. We conclude that mi/mg = 1/10 
is not within the scaling region and restrict the simultaneous fit to just data with mi/mg = 1/20. Fits to the = 12 
lattices for both the asqtad and HISQ/tree lattices are made at the single value of mi/mg = 1/20 where simulations 
were carried out. Our best fits to the iV^ = 12 asqtad data are shown in Fig. [19] for the 0(4) Ansatz. 

The final value of Tc is calculated by finding the location of the peak in the susceptibility. Errors on Tc for each 
action and Nt value are estimated by carrying out the entire analysis for 400 synthetic samples. Each point in these 
samples is taken from a Gaussian distribution with mean and standard deviation given by the central value and the 
quoted errors on the data points. Once the scaling function is determined, we can extract Tc at any desired value of 
mi/ms- Our final results at mi/mg = 1/20 and at the physical quark masses (defined as mi/mg = 0.037 with tuned 
to its physical value) are summarized in Table [Ij For the asqtad data, we extrapolate to to//tos = 0.0296 to partially 
correct for the fact that TOs on the LCP is about 20% heavier than the physical value. Changing the extrapolation 
point provides estimates at the correct physical light quark mass but does not correct for the heavier strange quark 
mass, so we expect the asqtad results to overestimate Tc- A final point, these estimates are also consistent with the 
location of the peak in xi,disc as shown in Figs. [TUl and [TT] and discussed in Sec. IIVBT] 

The last step in the determination of Tc is to extraplolate estimates at the physical value of mi/ms obtained at 
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different Nr to the continuum limit. With two asqtad data points and three HISQ/tree points, we explored all 
possible combinations of linear and quadratic fits to the data for each action and combined fits to the asqtad and 
HISQ/tree data in which the intercept (continuum value) is constrained to be the same for the two actions. The 
results of individual extrapolations linear in for the asqtad data and quadratic for the HISQ/tree data are 

given in rows marked oo in Table. HI The results of a combined quadratic fit are given in the final row marked "oo 
(asqtad+HISQ/tree)" . Two examples of combined fits are shown in Fig. [20l In Fig. [20l left) we show the result of a 
combined fit linear in the asqtad data and quadratic in the HISQ/trcc data, while the combined quadratic fit to both 
data sets is shown in Fig. EOfright). 

We take the result, T = 154(8) MeV, of the quadratic fit with the 0(4) Ansatz to the HISQ/tree data as our best 
estimate. These fits are less sensitive to variations in the range of temperature selected. The estimate from the 0(2) 
fits is about 2 MeV higher as shown in Table HI Since the quadratic fit to just the HISQ/tree data is identical to the 
fit shown in Fig. [201 right) we do not show it separately. For both 0(2) and 0(4) Ansatze, we find that the range of 
variation of the central value of the various fits is about 5 MeV, so 8 MeV is a conservative estimate of the combined 
statistical and systematic error and includes the full range of variation. Note that the error estimates, as expected, 
increase significantly for quadratic fits compared to linear fits, because in this case there are as many parameters as 
data points. In all these fits is much less than one. Furthermore, when the error bars on individual points do 
not represent normally distributed statistical errors and there are one or zero degrees of freedom, chi-squarc is not a 
useful guide for selecting the best fit. 

In determining this final estimate, we have mostly used the asqtad data as confirmatory for two reasons. First, the 
slope of the fits to the asqtad data is 2-2.5 times that in the HISQ/tree fits and the undetermined quadratic term 
may be large. Note that the large slope confirms our discussion in Sec. [H] that the discretization errors in the asqtad 
formulation are larger. Second, the LCP defined by the strange quark mass is ss 20% heavier than the physical value, 
so the asqtad action data overestimate Tc- 

In addition to the above error estimate, there is a 1 MeV uncertainty in the setting of the temperature scale as 
discussed in Sec. HI] which will shift all estimates. We therefore quote it as an independent error estimate. Thus, our 
result after continuum extrapolation, and at the physical light quark masses, mi/ms = 1/27, is Tc = (154±8± l)MeV. 
To obtain an overall single error estimate we add the uncertainty in the scale coming from ri to the statistical and 
systematic errors. Thus, our final estimate is 

Tc = (154 ± 9)MeV . (39) 



C. Comparison with previous results 

We compare the result in Eq. with three previous 2 + 1 flavor studies that also extrapolated Tc data to the 
continuum limit and to the physical light quark mass. It should be emphasized that in our previous work (20| such 
an extrapolation was not carried out. 
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Figure 17: Scaling fits and data for the chiral condensate Mb calculated with the HISQ/tree action on lattices with temporal 
extent A'^^ = 8 (left) and the chiral susceptibility Xm,!(right). The data for M(, at mi/nis = 0.025 and for and Xm,i at 
mi/nis = 0.05 are fit simultaneously using the 0(2) scaling Ansatz. The fits using the 0(4) Ansatz are similar. The points 
used in the scaling fits are plotted using open symbols. The dotted lines give the data scaled to the physical quark masses. 
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Figure 18: Scaling fits and data for tlie cliiral condensate Mb calculated with the HISQ/tree action on lattices with temporal 
extent Nt = 12 (left) and the chiral susceptibility Xm,; (right). The data for Mb and Xm.i at mi/nis — 0.05 are fit simultaneously 
using the 0(4) scaling Ansatz. The fits using the 0(2) Ansatz are similar. The points used in the scaling fits are plotted using 
open symbols. The dotted lines give the data scaled to the physical quark masses. 
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Figure 19: Scaling fits and data for the chiral condensate Mb calculated with the asqtad action on lattices with temporal extent 
Nt — 12 (left) and the chiral susceptibility Xm,; (right). The data for Mb and Xm.,i at mi/ms = 0.05 are fit simultaneously using 
the 0(4) scaling Ansatz. The fits using the 0{2) Ansatz are similar. The points used in the scaling fits are plotted using open 
symbols. The dotted lines give the data scaled to the physical quark masses. 



The 2004 study by the MILC collaboration used the asqtad action [T3| ■ They extrapolated Tc defined as the peak 
position in the total chiral susceptibility using an expression that incorporated the 0(4) critical exponent. They found 
Tc ~ 169(12)(4) in the chiral limit, which is just consistent, within errors, with our current result. Note, however, 
that the present data are more extensive and the scaling analysis is more comprehensive. 

The RBC-Bielefeld collaboration studied the p4 action on Nr — 4 and 6 lattices with several values of the light 



quark masses The result of a combined extrapolation in the quark mass and gave Tc = 192(7) (4) MeV 



This value is significantly higher than the one given in Eq. (j39p . Based on analyses done subsequently, and the new 
data on N^- = 8 lattices presented here, we find that this discrepancy is due to an underestimate of the slope of the 
linear fit, Tc vs. l/N^- The change in the slope between Nt = 4 and 6 data and Nt = 6 and 8 data is shown in 
Fig. 1211 To perform a continuum extrapolation of Tc obtained with the p4 action requires new calculations on lattices 
with Nt > 12. Having demonstrated that the discretization errors are similar for the p4 and asqtad actions and given 
the consistency between results obtained with the HISQ/tree, asqtad and stout actions, we do not intend to pursue 
further calculations with the p4 action. 

Lastly, the Wuppertal-Budapest collaboration has carried out a continuum extrapolation using the stout action 
on Nt = 6, 8, and 10 lattices, and more recently including Nt = 12 and 16 lattices (22| - [2^ . They work directly 
at mi/nis = 0.037 and perform a linear extrapolation in l/N"^ of their data to obtain continuum estimates. From 
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Figure 20: The data and the continuum extrapolation of the transition temperature obtained from the 0(4) scaUng analysis 
for mi /iris extrapolated to the physical value. The left figure shows a combined fit using a linear extrapolation to the asqtad 
data and a quadratic to the HISQ/tree data. The right figure shows a combined quadratic fit to both data sets. 
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180.6(4.7) 
167.4(4.2) 
156.8(8.4) 


181.9(3.2) 
170.1(3.1) 
160.7(6.1) 


184.5(4.4) 
172.5(3.9) 
162.9(7.9) 


185.0(3.4) 
174.0(2.9) 
165.2(5.9) 


6 (HISQ/tree) 

8 

12 

oo 


165.7(2.0) 
158.9(2.4) 
155.5(3.0) 
153.7(8.4) 


167.7(2.3) 
160.6(1.6) 
157.4(3.0) 
156.0(7.8) 


169.8(2.0) 
163.1(2.7) 
158.9(2.7) 
155.9(8.2) 


171.3(2.1) 
164.1(1.8) 
160.9(2.1) 
159.5(6.1) 


oo (asqtad + HISQ/tree) 


153.7(8) 


156.0(8) 


155.9(8) 


159.5(6) 



Table I: Pseudocritical temperature Tc determined from 0(2) and 0(4) scaling fits to the chiral condensate and the chiral 
susceptibility. The results are shown at mi/ms = 0.05 at which simulations have been performed for all three values of A''t and 
at the physical value of mi /iris. The temperature scale used in the fits is set using ri. Rows labeled oo give results after a 
linear extrapolation to the continuum for the asqtad data and quadratic for the HISQ data. The last row gives the results of 
a combined fit using an extrapolation quadratic in l/N'^, which coincides with the extrapolation of just the HISQ/tree data. 

the position of the peak in the renormalized chiral susceptibility they extract Tc = 147(2) (3) MeV. This value is 
approximately 1 a lower than our result. They also report higher transition temperatures of Tc — 157(3) (3) and 
155(3) (3) derived from inflection points in A/,s and {iIjiP)fi^ respectively. As discussed in Sec. Illll and confirmed by the 
significant difference in their three estimates, using inflection points is a less sensitive probe of the critical behavior, 
and thus less reliable for extracting Tc- By making simultaneous fits to the chiral condensate and the susceptibility, 
and defining the transition temperature as the peak in the susceptibility, our analysis overcomes this ambiguity. 

VI. QUARK NUMBER SUSCEPTIBILITIES AND THE POLYAKOV LOOP 

In this section, we present our results for light and strange quark number susceptibilities as well as the Polyakov 
loop expectation value. In Sec. IIIIl we stated that we do not use these observables for the determination of the QCD 
transition temperature because the singular part of the free energy, which drives the chiral transition, is subleading in 
the quark number susceptibilities and thus difficult to isolate. Similarly, in the case of the Polyakov loop expectation 
value, a direct relation to the critical behavior in the light quark mass regime has not been established. Nonetheless, 
we study these observables as they provide important insight into the deconfining aspects of the QCD transition. 

The discussion in Sec. IIIIl shows that quark number susceptibilities probe whether the relevant degrees of freedom 
of the system at a given temperature are hadronic or partonic; a rapid rise in the quark number susceptibilities 
signals the increasing contribution of light (partonic) degrees of freedom. Theoretically, the reduced temperature 
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Figure 21: Linear extrapolations to the continuum limit of Tc data for the p4 action. The fit to the Nt = 4 and 6 data is shown 
in red, while that to the Nt — 6 and 8 points is in blue. The two fits illustrate the change in the slope as Nt is increased. 
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Figure 22: Light quark number susceptibility for the HlSQ/tree action with mi/nis = 1/20 are compared with the strange 
quark number susceptibility. In the left panel ri is used to set the lattice scale, while in the right panel fx is used. The filled 
squares correspond to A^r = 12. 



t defined in Eq. (|10p is symmetric in the light and strange flavors at leading order, therefore, the structure of the 
singular contribution to both observables is similar and should differ only in magnitude. We compare the temperature 
dependence of xi and Xs on Nj- = 6 and 8 lattices in Fig. [52] and find a similar behavior with xi exhibiting a more 
rapid rise at temperatures closer to the chiral transition temperature. From these data, one would be led to conclude 
that xi E^iid Xs are probing the singular part of the free energy and deduce a higher transition temperature from Xs 
than from xi- However, in Sec. IIIII we pointed out that the temperature dependence of these observables, even in the 
chiral limit, is dominated by the regular part of the free energy. This feature is highlighted by the difference in the 
temperature dependence of xi and Xs versus that of the chiral condensate, shown in Fig. [51 which is dominated by the 
singular part in the chiral limit. Thus, even though the data in Fig. l22l show a rapid crossover, extracting information 
about the singular part from them is nontrivial, and consequently the determination of Tc from these observables is 
less reliable. 

In Fig. I22[ we also compare the temperature dependence of xi and Xs using ri and fx to set the temperature scale. 
Both light and strange quark number susceptibilities show a much smaller cutoff dependence when using the /k scale. 
A similar behavior was observed for the chiral susceptibility as discussed in Sec. IIV Bl 

In order to analyze further the cutoff dependence and to compare our data with the continuum extrapolated stout 
results [2^ . we examine the strange quark number susceptibility, which is statistically under better control. In Fig. 1231 
we show Xs/T"^ for m; = O.OStos calculated with the HISQ/trce and asqtad actions. The larger cutoff effects are in 
the asqtad data. We also plot the data using fx to set the temperature scale and find that all data collapse into 
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Figure 23: Data for the strange quark number susceptibility for mi = 0.05ms with the asqtad and HISQ/tree actions are 
compared to the continuum extrapolated stout results [SJ- In the right panel, we show that all the data collapse to a single 
curve when fx is used to set the temperature scale. 
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Figure 24: The renormalized Polyakov loop calculated for the asqtad and HISQ/tree actions and compared with the contirmum 
extrapolated stout result [S^]. In the right panel, we show the same data using /k to set the scale. 



one curve as shown in Fig. I23r right). suggesting that the discretization errors in Xs and /k are very similar. On the 
basis of this agreement, we conclude that the differences in results for Xs seen in Fig. I^STleft) are accounted for by 
discretization errors. 

The Polyakov loop is a sensitive probe of the thermal properties of the medium even though it is not a good 
measure of the critical behavior as discussed in Sec. IIIII After proper renormalization, the square of the Polyakov 
loop characterizes the long distance behavior of the static quark anti-quark free energy; it gives the excess free energy 
needed to screen two well separated color charges. The renormalized Polyakov loop Lren(T) has been studied in the 
pas t in the pure gauge theory [t^, [ill as well as in QCD with two [t^, three HO] and two plus one flavors of quarks 
|l9l[20l|. Following Ref. [22| Lren(T'), for a given Nr, is obtained from the bare Polyakov loop defined in Eq. <\22^ as 

Lren{T) = CXp{~ NrC{f3) / 2)Lbare , (40) 

where the renormalization constant c(/3) is deduced from the heavy quark potential discussed in Sec. HH 

We compare results for Lren{T) from the asqtad, HISQ/tree and stout actions in Fig. [M] to determine the cutoff 
effects in the description of basic thermal properties, e.g., in the asymptotic long distance structure of a heavy quark 
free energy and the screening of the quark- ant iquark force. We find significant cutoff effects in Lren with the asqtad 
action; however, these are considerably smaller than seen in the other quantities discussed so far. This is partially 
expected as the Polyakov loop is a purely gluonic observable and is thus less affected by the taste symmetry breaking 
in the staggered formulations. In Fig. [24l right), we also show the temperature dependence of Lren when fx is used 
to set the scale. For the HISQ/tree action, the already small cutoff effects are further reduced and the data move 
even closer to the continuum extrapolated stout result. For the asqtad data, on the other hand, the Nr dependence 
is slightly more pronounced. The more remarkable feature is the switch in the approach to the continuum limit; it is 
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Figure 25: The renormalized Polyakov loop for the HISQ/tree action on Nt — 6 lattices for two values of the light quark mass 
mt/m^ = 0.05 and 0.20. 



from above when /k is used to set the temperature scale and from below when using the ri scale. This is similar to 
what has been observed in calculations with the stout action on coarser lattices |22| . Such behavior is not surprising, 
as there is no obvious reason for cutoff effects in a hadronic observable like Jk and a gluonic observable like L^en to be 
related so that there is a cancellation of lattice artifacts in their ratio. Based on these observed changes with Nr and 
dependence on whether the scale is set by ri or fx, we expect that the asqtad action will need Lren data on Nr > 12 
lattices to agree with the continuum extrapolated HISQ/tree and stout data. This confirms the statement made in 
SeclUthat while asqtad and the HISQ actions are equivalent at order ©(a^a^) the difference at higher orders makes 
sizeable contributions below ^ 12. Overall, we conclude that once discretization effects are accounted for, the 
Polyakov loop measured using different staggered discretization schemes yields consistent estimates of the thermal 
properties of the theory. 

Lastly, we briefly discuss the quark mass dependence of the Polyakov loop. In Fig. [551 we show the renormalized 
Polyakov loop for two quark masses rrii/ms = 0.05 and 0.20 on Nt = 6 HISQ/tree ensembles. The quark mass 
dependence is small in the transition region, and at higher temperatures the renormalized Polyakov loop rises faster 
for the heavier quark mass. This dependence on the quark mass is consistent with the results obtained using the p4 
action [2lj. Previous work suggests that for larger values of the light quark masses, mi > 0.2ms, the quark mass 
dependence becomes more significant as L is influenced by the phase transition in the pure gauge theory [80|. 



VII. CONCLUSIONS 



In this paper, we present results for the chiral and deconfinement aspects of the QCD transition using three 
improved staggered quark actions called p4, asqtad and HISQ/tree. The chiral transition is studied using the chiral 
condensate and the chiral susceptibility, while the deconfining aspects of the transition have been addressed using the 
quark number susceptibilities and the renormalized Polyakov loop. By comparing the results obtained at different 
lattice spacings, we analyze the cutoff dependence of these quantities, in particular, the effects due to taste symmetry 
breaking. We find that the effects of taste symmetry breaking are significantly reduced in the HISQ/tree action. After 
demonstrating control over systematic errors, we compare our data with the continuum extrapolated stout results 
and find agreement for all quantities. 

We analyze the chiral transition in terms of universal 0{N) scaling functions using simulations performed at 
several values of the light quark mass and Nt- We show that for mi /iris < 0.1 the corrections to scaling are small and 
determine the critical temperature in the 2-fiavor (degenerate u and d quarks) chiral limit. We define the pseudocritical 
temperature Tc at non-zero values of the light quark mass, in particular the physical value mi/ms = 1/27, as the 
peak in the total chiral susceptibility. We find that this definition of the pseudocritical temperature gives estimates 
consistent with those obtained from the location of the peak in the disconnected chiral susceptibility used in previous 
studies. Our final result, after extrapolation to the continuum limit and for the physical value of light quarks, 
mi/mg = 1/27, is Tc = (154 ± 9)MeV based on the analysis summarized in Table ID 

The Wuppertal-Budapest collaboration has extracted Tc using the infiection point in the data for the chiral con- 
densate and the peak in the chiral susceptibility. In their most recent publication p^ . they quote values ranging 
from 147 to 157 McV depending on the chiral observable. While their methods for extracting Tc differ, their values 
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are, within errors, in agreement with our determination. The consistency of the results obtained with the HISQ/tree, 
asqtad and stout actions shows that these resuhs represent the continuum Hmit value for the staggered formulation. 

Finally, we also consider it interesting to determine the pseudocritical temperature in terms of other scaling func- 
tions, namely the scaling functions related to the mixed susceptibility, Eq. (HH]), and the specific heat. This would 
allow estimation of the intrinsic uncertainty in defining the pseudocritical temperature at nonzero values of the light 
quark mass where, a priori, it can depend on the probe used to estimate it. We will address this issue in future work. 
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mi 


= 0.05ms 


mi = O.lOma 


mi — 0.20ms 




nil 




2 X TUt 


mi 


ms 


2 X TUt 


mi 


ms 


2 X TUt 


3.4900 





001450 


0.0290 


31680 














3.5000 





001265 


0.0253 


32520 














3.5100 





001300 


0.0260 


30050 














3.5150 





001200 


0.0240 


30480 


— 


— 


— 


0.00480 


0.0240 


3510 


3.5200 





001200 


0.0240 


39990 


0.00240 


0.0240 


11800 


0.00480 


0.0240 


6470 


3.5225 





001200 


0.0240 


54610 


0.00240 


0.0240 


30620 


— 


— 


— 


3.5250 





001200 


0.0240 


64490 


0.00240 


0.0240 


31510 


— 


— 


— 


3.5275 





001200 


0.0240 


65100 


0.00240 


0.0240 


28060 


— 


— 


— 


3.5300 





001200 


0.0240 


70230 


0.00240 


0.0240 


34420 


0.00480 


0.0240 


24900 


3.5325 




— 


— 


— 


0.00240 


0.0240 


27310 


— 


— 


— 


3.5350 





001200 


0.0240 


39510 


0.00240 


0.0240 


29250 


0.00480 


0.0240 


4850 


3.5375 




— 


— 


— 


0.00240 


0.0240 


30320 


0.00480 


0.0240 


4890 


3.5400 





001200 


0.0240 


56740 


0.00240 


0.0240 


31590 


0.00480 


0.0240 


5190 


3.5425 




— 


— 


— 


0.00240 


0.0240 


38520 


0.00480 


0.0240 


4380 


3.5450 





001075 


0.0215 


20700 


0.00240 


0.0240 


15060 


0.00480 


0.0240 


5300 


3.5475 
















0.00480 


0.0240 


4630 


3.5500 
















0.00480 


0.0240 


5530 


3.5525 
















0.00480 


0.0240 


22540 


3.5550 
















0.00480 


0.0240 


23420 


3.5600 





001025 


0.0205 


11660 








0.00480 


0.0240 


21230 


3.5700 
















0.00480 


0.0240 


9250 



Table II: Parameters used in simulations with the p4 action on A^,- = 8 lattices and three lines of constant physics defined by 
mi/ms — 0.2, 0.1 and 0.05. The quark masses are given in units of the lattice spacing a. The finite temperature runs were 
carried out on 32^ x 8 lattices. The statistics are given in units of trajectories (2 x TUt) which are of length rum ~ 0.5. 

Appendix A: Lattice parameters for p4, asqtad and HISQ/tree simulations 

In this appendix, we provide details on the parameters used in the calculations with the p4, asqtad and HISQ/tree 
actions. These include the gauge couplings, quark masses and lattice sizes at which simulations were carried out and 
the collected statistics. In the following tables, TUt stands for the molecular dynamics time units simulated for finite 
temperature runs and TUq for the corresponding zero temperature ones. In Table [TTl we give the parameters of the 
simulations with p4 action on 32'^ x 8 lattices. In Table Hill and Table ITVl we give the simulation parameters for the 
asqtad action corresponding to Nr = 8 and Nr = 12 lattices, respectively. The parameters for the Nr — 8 ensemble 
with mi/ms = 0.1 runs are given in Ref. [20| . The details on the extraction of tq and ri are given in Sec. Ill Cl aud in 
Appendix iBl 

For the HISQ/tree action, the A^^- = 6 calculations with m; = 0.2ms were done on 16"^ x 32 and 16"^ x 6 lattices for 
zero and finite temperature studies, respectively, and the run parameters are summarized in Table [Vj 

The parameters of finite temperature simulations at m; = 0.05ms for ensembles with A^,- = 6, 8 and 12 and the 
corresponding zero temperature calculations are given in Table [VTl For the 0{N) scaling analysis, we have also used 
the HISQ/tree data on Nt = 6 and 8 lattices at mi = 0.025ms as summarized in Table |VlTl These data are part of 
the ongoing RBC-Biclcfcld study of 0{N) scaling [76j . 

Data for the light and strange quark chiral condensate, disconnected and connected chiral susceptibility, the light 
and strange quark number susceptibility and the bare Polyakov loop on finite temperature ensembles for the asqtad 
action on N^. — 8 and 12 lattices are given in Tables [VIIII and HXl The corresponding data for the HISQ/tree action 
are given m Tables 1x1 IXll and IXlIl for ^ 6, 8 and 12 lattices, respectively. 

For completeness, wc also give the data for the light and strange quark chiral condensate on zero temperature 
ensembles for the asqtad action on Nt = 8 and 12 lattices in Tabic IXlIIl The corresponding data for the HISQ/tree 
action arc given in Table IXIVI 
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mi - 


-- 0.05ms 




mi — 0.20m3 




nil 


rris 


Mo 


TUo 


TUt 


mi 


ma 




TUt 


6 


5500 


0.00353 


0.0705 


0.8594 


— 


29541 


— 




— 




— 


— 


6 


5750 


0.00338 


0.0677 


0.8605 


— 


13795 


0.01390 





0684 
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8515 


6 


6000 


0.00325 


0.0650 


0.8616 
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25979 


0.01330 
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8614 
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6 


6250 


0.00312 


0.0624 


0.8626 


3815 
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0.01270 





0628 





8624 


12695 


6 


6500 


0.00299 


0.0599 


0.8636 


3810 


22731 


0.01210 





0602 





8634 


12695 


6 


6580 


0.00295 


0.0590 


0.8640 


3789 


32725 


— 




— 




— 


— 


6 


6660 


0.00292 


0.0583 


0.8643 


3445 


32020 


— 




— 




— 


— 


6 


6750 


0.00288 


0.0575 


0.8646 


4840 


31067 


0.01160 





0577 





8645 


13695 


6 


6880 


0.00282 


0.0563 


0.8652 


3895 


33240 


— 










— 


6 


7000 


0.00276 


0.0552 


0.8656 


5475 


32324 


0.01110 





0552 





8655 


12695 


6 


7300 


0.00262 


0.0525 


0.8667 


3827 


26634 


0.01050 





0526 





8666 


12655 


6 


7600 


0.00250 


0.0500 


0.8678 


4725 


30433 


0.01000 





0500 





8677 


10695 


6 


8000 


0.00235 


0.0471 


0.8692 


5710 


35690 


0.00953 





0467 





8692 


7905 


6 


8300 












0.00920 





0452 





8702 


6910 



Table III: Parameters used in simulations with the asqtad action on A'^^ = 8 lattices and two lines of constant physics defined 
by mi /ma = 0.05 and 0.2. The quark masses are given in units of the lattice spacing a. The zero and finite temperature runs 
were carried out on 32* and 32'^ x 8 lattices, respectively. The number of molecular dynamics time units for zero and finite 
temperature runs are given in columns labeled by TUo and TUt, respectively. 





mi 


ma 


Mo 


TUo 


TUt 


6 


800 


0.00236 


0.0471 





8692 


2040 


29687 


6 


850 


0.00218 


0.0436 





8709 


3745 


30452 
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875 


0.00210 


0.0420 





8718 


2025 


18814 
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900 


0.00202 


0.0404 





8726 


3340 


33427 
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925 


0.00195 


0.0389 





8733 


2925 


27966 
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950 


0.00187 


0.0375 





8741 


2890 


27734 
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975 


0.00180 


0.0361 





8749 




28464 
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985 


0.00178 


0.0355 
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25191 
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000 


0.00174 


0.0347 





8756 


3055 


25580 
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025 


0.00167 


0.0334 





8764 




21317 
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050 


0.00161 


0.0322 
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085 


0.00153 


0.0305 





8781 
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0.00144 


0.0289 





8793 


3100 


32640 


7 


200 


0.00128 


0.0256 





8813 


4615 


31175 


7 


300 


0.00110 


0.0220 





8839 


2645 


34240 


7 


400 


0.000946 


0.0189 





8863 


1611 


22805 


7 


550 


0.000754 


0.0151 





8863 


1611 


26415 



Table IV: Parameters used in simulations with the asqtad action on Nt = 12 lattices and the LCP defined by mi /ma = 0.05. 
The quark masses are given in units of the lattice spacing a. The zero and finite temperature runs were carried out on 48'* and 
48^ X 12 lattices respectively. 



Appendix B: Setting the lattice spacing in asqtad and HISQ/tree simulations 

Three obscrvablcs, tq, ri and fx, were studied to set the lattice scale as discussed in Sec. [Ill In simulations with the 
asqtad action, we used the values of ri published in Ref. [4l| and performed an interpolation using a renormalization 
group inspired Ansatz j4ll | 

a _ co/(/3)+C2(10//3)/3(/3)+e4(10//3)V3(/3) 

n l + rf2(10//?)/2(/3)+d4(10//3)2/4(/?) ' ^ > 
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TUt 
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000 





0230 





115 


3,000 


6,000 


6 


038 





0216 





108 


3,000 


6,000 


6 


100 





0200 





100 


3,000 


6,000 


6 


167 





0182 





091 


3,000 


6,000 


6 


200 





0174 





087 


3,000 


6,000 


6 


227 





0168 





084 


3,000 


6,000 


6 


256 





0162 





081 


3,000 


6,000 


6 


285 





0158 





079 


3,000 


6,000 


6 


313 





0152 





076 


3,000 


6,000 


6 


341 





0148 





074 


3,000 


6,000 


6 


369 





0144 





072 


3,000 


6,000 


6 


396 





0140 





070 


3,000 


6,000 


6 


450 





0136 





068 


3,000 


6,000 


6 


800 





0100 





050 


3,000 





Table V: Parameters used in simulations with the HISQ/tree action on Nt = 6 lattices and the LCP defined by mi/nis = 0.2. 
The quark masses are given in units of the lattice spacing a. The zero and finite temperature runs were carried out on 16-' X 32 
and 16'^ x 6 lattices respectively. 



/(/3) = (6o(10//?))-''^/(2''°) exp(-/3/(206o)), (B2) 

Co = Coo + {coium + coisTOs + co2(2to/ + rris))/ J{I3) (B3) 

C2 = C2o + C2i(2mj + m,,)//(^). (B4) 

Here 6o and hi are the coefficients of the universal 2-loop beta functions and for the numerical values of the other 
coefficients we get 

Coo = 4.574615 • 10\ coi„ = 6.081198 • lO^i, cqi^ = 2.689340 • lO^^, 
co2 = -3.591183 • 10-3, C20 = -5.36 8 781 • 10^ C21 = 8.756186 • 10^, 

C4 = 2.930565 • 10^ ^2 = -3.786570 • 10^, ^4 = 7.385881 • 10^ (B5) 

To set the lattice scale using one needs to measure qJk on lattices with large volumes and on high statistics 
ensembles. The zero temperature asqtad ensembles generated in this study (see Tables HII] and HV]) were primarily 
for performing subtractions of UV divergences. We, therefore, used the MILC collaboration study of the light meson 
decay constants and masses [8l|. A systematic fit of fx data, obtained on a large set of ensembles with similar 
parameters values as in this study, using staggered chiral perturbation theory gives values for fxri at nine j3 values 
in the range [6.458,7.47]. We interpolated these results to obtain Jk^i at the physical quark mass at the /3 values 
used in the 0.0577is LCP study and given in Tables Hill and HVl With fK^i and our parametrization of ri/a in hand 
we calculated the temperature associated with each (3. The uncertainty in these T values is less than 1% over the 
range of j3 values investigated. 

In simulations with the HISQ/tree action, we used tq, ri and Jk to set the lattice spacing. We determined ro/a 
and ri/a values from the static quark anti-quark potential. These values are summarized in Tables IXVI and IXVII 
for the m; = 0.2ms and the m; = 0.05ms LCP, respectively. The last column contains the values of the additive 
renormalization constant c(/3) defined in Eq. and used in the calculation of the renormalized Polyakov loop. 

For the m; = 0.2ms LCP, we ffi r^/a values with a form similar to Eq. (|Bip : 

^10. _ cof{(3) +02(10/ (3) f (13) 

^^l/^V,=o.2™. - i + d2{W/f3)p{(3) ' 

where cq = 32.83, C2 = 81127, d2 = 1778 and xV^of = l-Ol- Then we use tq = 0.468 fm obtained in Sec. [Tl] to 
calculate the lattice spacing in units of fm. 

For the 0.05ms ensembles, we preferred to use the ri scale, since the value of ri in fm is accurately determined by 
the MILC collaboration [H^ and serves as an external input for this work. As one can see from Table IXVTl our data 
set includes coarse lattices where ri/a < 2. Extracting the value of ri/a at such short distances is problematic due to 
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Table VI: Parameters used in simulations with the HISQ/tree action on Nr = 6, 8 and 12 lattices and the LCP defined by 
mi/ms = 0.05. The quark masses are given in units of the lattice spacing a. The statistics in molecular dynamics time units 
TU are given for both the zero and finite temperature runs. The lattice sizes used for the N^ = 6, 8 and 12 finite temperature 
simulations were 24=^ X 6, 32=* X 8 and 48=* x 12, respectively. 
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Table VII: Parameters used in simulations with the HISQ/tree action on A'^ = 6 and 8 lattices and the LCP defined by 
mi I ma — 0.025. The quark masses are given in units of the lattice spacing a. The finite temperature runs were carried out on 
32^ X 6 and 32^ x 8 lattices as part of the RBC-Bielefeld study of 0{N) scaUng 



significant systematic errors. Therefore, we used the following strategy: for coarser lattices we converted rg/a values 
into ri/a using the continuum ratio (ro/ri)cont derived in Sec. HIl and defined as 



a 



^ \ ro/a/{ro/ri)cont, /? < /3oi, (-gy^ 

mi=0.05m, 1 Ti /fl, /3 ^ /3oi , 



where the optimal value of /3oi is determined as follows. We fit a/ri to the Ansatz: 

and varied /3oi in the range [6.423, 6.608] examining the of the fit. The minimum of is achieved at /3oi = 6.423. 
With this choice of /3oi the coefficients of the fit are cq = 44.06, ca = 272102, da = 4281 and x^/dof = 0.31. The data 
and the fit are shown in Fig. [26l 

We also use fx as an independent observable for scale setting. We fit data to a form similar to (|B8[) : 



In the continuum limit, the product rifx is fixed: 



0.3106 fm • 156.1/V2 MeV , 
197.3 fm.MeV ^ ^^^^^ 



Taking the ratio of Eqs. fB9| and (|B8| one finds 



ri.fK = ^ => c^ = 7M. (Bll) 
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Table VIII: Data for the light and strange quark chiral condensate, disconnected and connected chiral susceptibility, the strange 
quark number suceptibihty and the bare Polyakov loop for the asqtad action on A^^ = 8 lattices and the two lines of constant 
physics defined by mi/ms = 0.05 and mi/rris = 0.1. The data for the light and strange quark condensate are presented using 
the 2-flavor normahzation. 



We, therefore, varied only and df^ in the Ansatz (|B9p and got cf^ = 32911, = 2388, x^/dof = 3.5. Presumably, 
the errors on afx are somewhat underestimated and this results in a high x^/dof. However, the fluctuations of the 
data around the fit are, at worst, 1%; therefore we feel comfortable using the Ansatz in Eq. (|B9p for setting the 
temperature scale. This fit is shown in Fig. [27l 

For those ensembles where a/ri and u/k measurements are available, we can compare the difference between the 
lattice spacing defined from ri and fx- We define 
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Table IX: Data for the light and strange quark chiral condensate, disconnected and connected chiral susceptibility, the strange 
quark number susceptibility and the bare Polyakov loop for the asqtad action on Nt = 12 lattices and the LCP defined by 
nii/ms = 0.05. The data for the light and strange quark condensate are presented using the 2-flavor normalization. 
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Figure 26: The a/ri data, defined in Eq. (|B7|) . together with the smoothing fit for the HISQ/trec action, m; — 0.05ms LCP. 

and 

A, = 1°^" • 100%. (B14) 

an 

The difference Aa is shown in Fig.[28l It is about 8% in the lattice spacings corresponding to the coarsest ensemble and 
decreases to about 1% for the finest N^- = 12 ensembles on which fx was measured. The temperature corresponding 
to these fine lattices is about 240 MeV. 



Appendix C: Hadron spectrum with the HISQ/tree action 

This appendix contains results of the calculations of the hadron correlation functions and some details on the 
extraction of the masses and decay constants. Fits to the correlation functions for mesons and baryons are performed 
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0.134460(098) 


1.3652(617) 


1.2792(091) 


0.531(17) 


0.2740(40) 


0.02875(36) 


6.175 


174.610 


0.0906 


0.024075(159) 


0.127619(079) 


1.1297(692) 


1.3365(081) 


0.600(12) 


0.3200(54) 




6.195 


178.054 


0.0880 


0.020080(104) 


0.122064(051) 


1.0001(172) 


1.3415(068) 


0.640(35) 


0.3560(72) 


0.03664(29) 


6.215 


181.568 


0.0862 


0.016410(172) 


0.117038(095) 


0.7007(439) 


1.3053(074) 


0.675(11) 


0.3974(58) 


0.04057(31) 


6.245 


186.969 


0.0830 


0.012735(179) 


0.109962(123) 


0.4768(400) 


1.2076(075) 


0.711(07) 


0.4424(40) 


0.04646(49) 


6.285 


194.423 


0.0790 


0.009099(091) 


0.100599(080) 


0.1752(094) 


1.0271(050) 


0.757(14) 


0.5116(65) 


0.05534(40) 


6.341 


205.355 


0.0740 


0.006545(101) 


0.089489(156) 


0.0812(104) 


0.8432(088) 


0.799(10) 


0.5879(72) 


0.06806(54) 


6.354 


207.977 


0.0728 


0.006011(038) 


0.086944(096) 


0.0476(049) 


0.7958(056) 


0.798(10) 


0.6098(72) 


0.07110(54) 


6.423 


222.450 


0.0670 


0.004457(023) 


0.075293(079) 


0.0174(036) 


0.6583(028) 


0.838(08) 


0.6883(54) 


0.08810(38) 


6.488 


236.960 


0.0620 


0.003643(019) 


0.066299(040) 


0.0065(022) 


0.5854(017) 


0.846(06) 


0.7420(43) 


0.10367(41) 


6.515 


243.246 


0.0603 


0.003387(013) 


0.063337(051) 


0.0019(012) 


0.5632(016) 


0.849(05) 


0.7560(43) 


0.11238(52) 


6.550 


251.627 


0.0581 


0.003149(019) 


0.059679(031) 


0.0003(008) 


0.5410(012) 


0.849(08) 


0.7787(54) 


0.12180(61) 


6.575 


257.777 


0.0564 


0.002977(011) 


0.057089(046) 


0.0019(010) 








0.12714(58) 


6.608 


266.106 


0.0542 


0.002798(011) 


0.053897(029) 


0.0021(015) 








0.13663(63) 


6.664 


280.807 


0.0514 


0.002549(017) 


0.049855(046) 


0.0019(011) 


0.4966(009) 






0.15203(70) 


6.800 


319.609 


0.0448 


0.002101(012) 


0.041501(018) 


0.0005(007) 


0.4669(003) 


0.861(04) 


0.8248(36) 


0.18747(68) 


6.950 


367.875 


0.0386 


0.001746(007) 


0.034571(012) 


0.0004(009) 


0.4495(003) 


0.864(06) 


0.8374(58) 


0.22876(74) 


7.150 


442.160 


0.0320 


0.001412(010) 


0.027810(016) 


-0.0025(014) 


0.4353(001) 


0.849(05) 


0.8395(43) 


0.28241(73) 



Table X: Data for the light and strange quark chiral condensate, disconnected and connected chiral susceptibility, the light and 
strange quark number susceptibility and the bare Polyakov loop for the HISQ/tree action on A^,- — 6 lattices and the two lines 
of constant physics defined by mi/ms = 0.025 and mi /iris = 0.05. A thousand trajectories were discarded for thermalization. 

choosing either the maximum allowed separation, typically Nt/2 = 16, or up to a distance where the fractional error 
exceeds 30%. The minimum distance is varied to estimate the systematic effects due to contamination from excited 
states. In most cases, the fit function includes the ground state and the opposite parity state that contributes with 
an alternating sign, and their contribution due to back-propagation. In some cases, the first excited state is also 
included. To reduce the uncertainty due to auto-correlations in the estimation of statistical errors, we block the data 
by 10. 

The calculation of hadron correlators was done using wall sources. The pseudoscalar meson masses were calculated 
on most of our zero temperature lattices in order to establish the LCP. The data for the m; = 0.2m., LCP are given in 
Table IXVlIl and for O.OSm,, in Table IXVlIII We also calculated the pseudoscalar meson decay constants for different 
combinations of the quark masses. For these calculations, we used both wall and point sources and performed a 
simultaneous fit to the two correlators. In most cases, we used one wall source per lattice, but in some cases we used 
two or four sources. In Table IXIXl we give the pseudoscalar meson decay constants in lattice units and specify the 
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T MeV 


anis 


2 X (V-V)! 


2 X {iIjiIj}s 


dis 

Xi 


con 

Xi 


Xi 


Xs 


3 X -/vbare 


mi /iris = 0.025 


6.2450 
6.2600 
6.2850 
6.3000 
6.3150 
6.3300 
6.3540 
6.3650 
6.3900 
6.4230 
6.4450 


140.227 
142.298 
145.817 
147.970 
150.154 
152.370 
155.983 
157.666 
161.558 
166.837 
170.448 


0.0830 
0.0811 
0.0790 
0.0773 
0.0759 
0.0746 
0.0728 
0.0717 
0.0695 
0.0670 
0.0653 


0.03034(07) 
0.02812(14) 
0.02514(15) 
0.02279(19) 
0.02105(15) 
0.01953(17) 
0.01610(20) 
0.01427(17) 
0.01142(26) 
0.00842(20) 
0.00646(25) 


0.11810(4) 
0.11447(7) 
0.10994(8) 
0.10643(11) 
0.10381(7) 
0.10111(8) 
0.09688(8) 
0.09449(6) 
0.09010(9) 
0.08503(9) 
0.08143(15) 


0.90(07) 
1.09(10) 
0.96(09) 
1.16(11) 
1.30(07) 
1.32(11) 
1.47(07) 
1.44(07) 
1.32(09) 
0.77(10) 
0.72(10) 










rni/ms — 0.05 


6.1950 
6.2450 
6.2600 
6.2850 
6.3150 
6.3410 
6.3540 
6.3900 
6.4230 
6.4450 
6.4600 
6.4880 
6.5150 
6.5500 
6.5750 
6.6080 
6.6640 
6.8000 
6.9500 
7.1500 


133.541 
140.227 
142.298 
145.817 
150.154 
154.016 
155.983 
161.558 
166.837 
170.448 
172.952 
177.720 
182.435 
188.720 
193.333 
199.580 
210.605 
239.706 
275.906 
331.620 


0.0880 
0.0830 
0.0811 
0.0790 
0.0759 
0.0740 
0.0728 
0.0695 
0.0670 
0.0653 
0.0642 
0.0620 
0.0603 
0.0582 
0.0564 
0.0542 
0.0514 
0.0448 
0.0386 
0.0320 


0.042262(031) 
0.035540(065) 
0.033581(085) 
0.030407(073) 
0.026587(105) 
0.023661(103) 
0.021917(088) 
0.017638(160) 
0.013908(119) 
0.011433(130) 
0.010222(097) 
0.008284(082) 
0.006738(062) 
0.005356(051) 
0.004616(022) 
0.003944(024) 
0.003242(020) 
0.002343(009) 
0.001855(005) 
0.001454(003) 


0.12903(2) 
0.11837(4) 
0.11487(4) 
0.11018(3) 
0.10415(4) 
0.09968(4) 
0.09721(5) 
0.09054(5) 
0.08525(5) 
0.08153(7) 
0.07919(5) 
0.07519(5) 
0.07165(5) 
0.06722(5) 
0.06400(3) 
0.06007(4) 
0.05494(5) 
0.04462(2) 
0.03651(1) 
0.02901(1) 


0.3855(121) 
0.5719(183) 
0.5822(301) 
0.6314(262) 
0.7213(276) 
0.7739(205) 
0.7761(500) 
0.9397(436) 
0.8828(466) 
0.7199(464) 
0.6721(343) 
0.4976(242) 
0.3083(153) 
0.1661(165) 
0.0994(078) 
0.0636(043) 
0.0317(027) 
0.0068(011) 
-0.0013(007) 
0.0001(009) 


0.7987(25) 
0.8203(52) 
0.8395(42) 
0.8638(44) 
0.8978(56) 
0.9205(58) 
0.9320(53) 
0.9839(55) 
1.0351(67) 
1.0390(49) 
1.0315(65) 
0.9822(39) 
0.9126(44) 
0.8202(54) 
0.7540(65) 
0.6891(47) 
0.6125(26) 
0.5192(13) 
0.4786(06) 
0.4554(04) 


0.1670(314) 
0.2016(623) 
0.2559(183) 
0.2729(078) 
0.3481(175) 
0.3418(370) 
0.4147(230) 
0.4721(104) 
0.5547(051) 
0.6121(157) 
0.6473(043) 
0.6914(074) 
0.7296(100) 
0.7632(037) 
0.7878(114) 
0.8030(036) 
0.8318(031) 
0.8741(021) 
0.8984(016) 
0.9123(014) 


0.0659(51) 
0.0935(61) 
0.1095(28) 
0.1212(19) 
0.1461(50) 
0.1651(90) 
0.1849(64) 
0.2381(41) 
0.2895(24) 
0.3362(29) 
0.3623(30) 
0.4134(59) 
0.4589(58) 
0.5177(31) 
0.5632(64) 
0.6052(32) 
0.6711(29) 
0.7796(21) 
0.8443(17) 
0.8829(14) 


0.00289(06) 
0.00381(07) 
0.00406(10) 
0.00489(09) 
0.00577(06) 
0.00664(08) 
0.00734(06) 
0.00937(18) 
0.01151(10) 
0.01325(17) 
0.01468(11) 
0.01631(36) 
0.01937(16) 
0.02264(38) 
0.02539(20) 
0.02907(30) 
0.03566(28) 
0.05368(24) 
0.07535(24) 
0.10485(30) 



Table XI: Data for the light and strange quark chiral condensate, disconnected and connected chiral susceptibility, the light 
and strange quark number susceptibility and the bare Polyakov loop for the HISQ/tree action on A^^ ~ 8 lattices and the 
two lines of constant physics defined by mi/ms — 0.025 and mi/ms — 0.05. Fifteen hundred trajectories were discarded for 
thermalization. 



number of sources per lattice used in the calculations. 

We also calculated the vector meson and the baryon (nucleon and il) masses on some of the T = ensembles. The 
results are given in Table IXXl for the m; = 0.2ms LCP and in Table IXXJ for the mi — 0.05?tis LCP and discussed 
in Sec. |TT1 To improve the statistics we used 4 sources per lattice for the calculation of vector mesons and nucleon 
correlation functions. It was sufficient to use one source for the baryon as it is composed of three heavier strange 
quarks. 



Appendix D: Autocorrelations in HISQ/tree simulations 

An analysis of autocorrelations in zero-temperature calculations with the HISQ action at the light quark mass 
mi = 0.2ms has beeii presented in Ref. [i^. Note that the HISQ/tree action used in this paper differs from the 
HISQ action in Ref. [46| in that it docs not include a dynamical charm quark and uses only the tree-level instead 
of the one-loop improved coefficient for the Symanzik gauge action. We calculate the dimensionless autocorrelation 



38 



/5 


r MeV 


anis 


2 X (^V)i 


2 X 


dis 
Xl 


con 

Xl 


Xl 


Xs 


3 X Lhare 


6.6640 


140 


403 


0.0514 


0.011032(61) 


0.060172(32) 


0.414(57) 


0.6697(075) 


0.251(57) 


0.1440(072) 


0.001011(095) 


6.7000 


145 


321 


0.0496 


0.009642(38) 


0.057132(18) 


0.437(34) 


0.6652(155) 


0.310(52) 


0.1570(086) 


0.001173(110) 


6.7400 


150 


967 


0.0476 


0.008127(41) 


0.053835(28) 


0.449(23) 


0.6735(060) 


0.403(41) 


0.1987(115) 


0.001434(060) 


6.7700 


155 


329 


0.0461 


0.006998(78) 


0.051354(19) 


0.504(42) 


0.6748(110) 


0.504(37) 


0.2362(072) 


0.001814(075) 


6.8000 


159 


804 


0.0448 


0.006157(63) 


0.049359(18) 


0.543(51) 


0.6731(095) 


0.520(23) 


0.2606(072) 


0.002152(030) 


6.8400 


165 


949 


0.0430 


0.004685(76) 


0.046456(35) 


0.491(27) 


0.6921(068) 


0.657(24) 


0.3571(072) 


0.002821(101) 


6.8800 


172 


302 


0.0412 


0.003745(59) 


0.043759(28) 


0.391(39) 


0.6592(065) 


0.649(19) 


0.3816(101) 


0.003449(067) 


6.9100 


177 


207 


0.0401 


0.003099(47) 


0.041876(36) 


0.253(30) 


0.6421(121) 


0.757(14) 


0.4752(072) 


0.004293(106) 


6.9500 


183 


938 


0.0386 


0.002680(32) 


0.039783(26) 


0.164(13) 


0.6150(067) 


0.743(21) 


0.5069(086) 


0.005024(075) 


6.9900 


190 


893 


0.0370 


0.002307(36) 


0.037522(28) 


0.104(16) 


0.5743(034) 


0.814(17) 


0.5832(086) 


0.006139(092) 


7.0300 


198 


078 


0.0356 


0.002036(14) 


0.035543(24) 


0.062(07) 


0.5318(046) 


0.842(12) 


0.6322(101) 


0.007371(036) 


7.1000 


211 


225 


0.0333 


0.001738(10) 


0.032396(16) 


0.030(04) 


0.5052(046) 


0.855(09) 


0.6998(101) 


0.009511(077) 


7.1500 


221 


080 


0.0320 


0.001611(08) 


0.030768(08) 


0.016(03) 


0.4955(017) 


0.868(05) 


0.7387(043) 


0.011144(095) 


7.2800 


248 


626 


0.0285 


0.001361(07) 


0.026477(08) 


0.012(03) 


0.4710(007) 


0.891(04) 


0.8064(043) 


0.016284(077) 



Table XII: Data for the light and strange quark chiral condensate, disconnected and connected chiral susceptibility, the hght 
and strange quark number susceptibihty and the bare Polyakov loop for the HISQ/tree action on Nt — 12 lattices and the 
LCP defined by mi/nis = 0.05. Two hundred trajectories were discarded for thermalization. 



/3 




2 X (V-V); 


2 X 


32^ 


6.600 


0.0650 


0.04102(5) 


0.13908(4) 


6.625 


0.0624 


0.03716(6) 


0.13127(5) 


6.650 


0.0599 


0.03389(5) 


0.12406(4) 


6.658 


0.0590 


0.03266(7) 


0.12147(5) 


6.666 


0.0583 


0.03179(6) 


0.11944(3) 


6.675 


0.0575 


0.03086(4) 


0.11717(3) 


6.688 


0.0563 


0.02936(4) 


0.11375(3) 


6.700 


0.0552 


0.02811(5) 


0.11071(3) 


6.730 


0.0525 


0.02506(7) 


0.10331(4) 


6.760 


0.0500 


0.02260(5) 


0.09672(3) 


6.800 


0.0471 


0.01976(5) 


0.08896(2) 


48^ 


6.800 


0.0471 


0.01970(3) 


0.08886(2) 


6.850 


0.0436 


0.01672(4) 


0.07997(3) 


6.875 


0.0420 


0.01539(3) 


0.07594(2) 


6.900 


0.0404 


0.01425(2) 


0.07213(2) 


6.925 


0.0389 


0.01319(3) 


0.06855(2) 


6.950 


0.0375 


0.01218(2) 


0.06521(2) 


7.000 


0.0347 


0.01041(2) 


0.05888(1) 


7.050 


0.0322 


0.00898(3) 


0.05339(2) 


7.085 


0.0305 


0.00815(3) 


0.04982(2) 


7.125 


0.0289 


0.00730(2) 


0.04640(1) 


7.200 


0.0256 


0.00583(2) 


0.03989(1) 


7.300 


0.0220 


0.00447(3) 


0.03307(2) 


7.400 


0.0189 


0.00341(3) 


0.02748(1) 


7.550 


0.0151 


0.00223(5) 


0.02100(2) 



Table XIII: Data for the light and strange quark chiral condensate for the asqtad action on Nt = 8 and 12 zero temperature 
ensembles with the line of constant physics defined by mi/nia = 0.05. The data for the light and strange quark condensate are 
presented using the 2-flavor normalization. 
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ma 


2 X (V"^), 


2 X (V'V')s 


24^ X 32 


5.900 


0.1320 


0.09833(6) 


0.21164(3) 


6.000 


0.1138 


0.07801(22) 


0.18024(10) 


6.050 


0.1064 


0.06890(8) 


0.16621(6) 


28^ X 32 


6.100 


0.0998 


0.06070(7) 


0.15300(6) 


32-* 


6.195 


0.0880 


0.04700(14) 


0.12983(8) 


6.285 


0.0790 


0.03660(11) 


0.11134(3) 


6.341 


0.0740 


0.03128(8) 


0.10134(4) 


6.354 


0.0728 


0.03014(8) 


0.09900(5) 


6.390 


0.0695 


0.02723(7) 


0.09295(4) 


6.423 


0.0670 


0.02499(7) 


0.08817(3) 


6.488 


0.0620 


0.02069(9) 


0.07906(6) 


6.515 


0.0603 


0.01920(10) 


0.07588(5) 


6.550 


0.0582 


0.01744(9) 


0.07203(4) 


6.575 


0.0564 


0.01635(4) 


0.06908(3) 


6.608 


0.0542 


0.01485(6) 


0.06548(5) 


6.664 


0.0514 


0.01281(6) 


0.06061(3) 


6.800 


0.0448 


0.00883(6) 


0.05008(4) 


6.950 


0.0386 


0.00611(9) 


0.04130(3) 


7.150 


0.0320 


0.00366(8) 


0.03259(5) 


32^ X 64 


6.460 


0.0642 


0.02238(5) 


0.08269(2) 


48^ 


6.740 


0.0476 


0.01080(4) 


0.05445(2) 


6.880 


0.0412 


0.00765(3) 


0.04500(1) 


7.030 


0.0356 


0.00544(4) 


0.03731(1) 


48^ X 64 


7.150 


0.0320 


0.00420(5) 


0.03255(2) 


7.280 


0.0285 


0.00322(2) 


0.02815(1) 



Table XIV: Data for the light and strange quark chiral condensate for the HISQ/tree action on = 6, 8 and 12 zero 
temperature ensembles with the line of constant physics defined by mi/ras — 0.05. The data for the light and strange quark 
condensate are presented using the 2-flavor normalization. 



coefRcient defined as 



c^, ^ i-^-^+^t) - {x.y ^^^^ 



for several representative ensembles both at zero and finite temperature at m; = 0.05ms. 

The autocorrelation coefficients for the plaquctte, fight and strange quark chiral condensate and the topological 
charge are shown in Tables IXXIII and IXXIIII At zero temperature, the autocorrelation coefficient is given at time 
separation At = 5 for /3 = 6.000, 6.285 and 6.460, and At = 6 for /3 = 7.280. At finite temperature, the autocorrelation 
coefficient is given at time separation At = 10. The chiral condensate was measured using one source on zero 
temperature lattices and with ten sources on finite temperature lattices. In several cases (mostly for zero temperature 
ensembles), the autocorrelation function for the light quark chiral condensate is too noisy to extract Ca* reliably. 

The integrated autocorrelation time, 

oo 

T„,t = 1 + 2 ^ C(At) , (D2) 

At=l 



provides an estimate of the time separation at which measurements can be considered statistically independent. To 
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p 


ro/a 


ri/a 


fit 1 


c(/3) ■ ro 


6.000 


2.037(12) 


1.410(13) 


2.052 


-1.622(17) 


6.038 


2.141(12) 


1.473(12) 


2.128 


-1.678(28) 


6.100 


2.250(12) 


1.544(19) 


2.256 


-1.834(44) 


6.167 


2.413(12) 


1.659(19) 


2.403 


-1.965(66) 


6.200 


2.501(17) 


1.722(20) 


2.478 


-2.008(26) 


6.227 


2.537(12) 


1.745(12) 


2.542 


-2.094(52) 


6.256 


2.603(14) 


1.798(20) 


2.611 


-2.124(53) 


6.285 


2.715(25) 


1.813(60) 


2.683 


-2.167(32) 


6.313 


2.742(20) 


1.848(12) 


2.753 


-2.274(39) 


6.341 


2.802(32) 


1.910(10) 


2.826 


-2.299(54) 


6.369 


2.916(39) 


1.983(14) 


2.900 


-2.331(66) 


6.396 


2.937(30) 


2.016(20) 


2.973 


-2.420(61) 


6.450 


3.110(30) 


2.132(14) 


3.124 


-2.546(84) 


6.800 


4.330(65) 


2.962(26) 


4.278 


-3.458(73) 



Table XV: Estimates of the scale setting parameters ro, ri and the additive renormalization constant c(/3) in the determination 
of the potential for the HISQ/tree calculations along the m; = 0.2ms LCP. See text for the definition of rjj**. 



/3 


ro/a 


n/a 


fit 1 

ri fa 


c(/3) ■ ro 


5.900 


1.909(11) 


1.23(13) 


1.263 


-1.441(15) 


6.000 


2.094(21) 


1.386(80) 


1.391 


-1.639(28) 


6.050 


2.194(22) 


1.440(31) 


1.460 


-1.748(28) 


6.100 


2.289(21) 


1.522(30) 


1.533 


-1.828(27) 


6.195 


2.531(24) 


1.670(30) 


1.682 


-2.072(31) 


6.285 


2.750(30) 


1.822(30) 


1.836 


-2.257(36) 


6.341 


2.939(11) 


1.935(30) 


1.940 


-2.440(14) 


6.354 


2.986(41) 


1.959(30) 


1.964 


-2.498(49) 


6.423 


3.189(22) 


2.096(21) 


2.101 


-2.653(27) 


6.460 


3.282(32) 


2.165(20) 


2.178 


-2.706(36) 


6.488 


3.395(31) 


2.235(21) 


2.238 


-2.808(37) 


6.550 


3.585(14) 


2.369(21) 


2.377 


-2.946(17) 


6.608 


3.774(20) 


2.518(21) 


2.513 


-3.070(27) 


6.664 


3.994(14) 


2.644(23) 


2.652 


-3.251(16) 


6.800 


4.541(30) 


3.025(22) 


3.019 


-3.675(31) 


6.880 


4.901(18) 


3.246(22) 


3.255 


-3.896(18) 


6.950 


5.249(20) 


3.478(23) 


3.475 


-4.077(40) 


7.030 


5.668(49) 


3.728(26) 


3.742 


-4.439(47) 


7.150 


6.275(39) 


4.177(31) 


4.176 


-4.791(37) 


7.280 


6.991(72) 


4.705(26) 


4.697 


-5.210(89) 



Table XVI: Estimates of the scale setting parameters ro, ri and the additive renormalization constant c(/3) in the determination 
of the potential for the HISQ/tree calculations along the mi = 0.05ms LCP. See the text for the definition of r{^^ ja and fit 
details. 



calculate Tint reliably requires time series substantially longer than generated in this study, so we provide a rough 
estimate. We assume that the autocorrelation function is dominated by a single exponential, in which case 



C{M) = exp(-AVTi) 



(D3) 



and 



Tl 



At 



(D4) 
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6.000 


0.36320(15) 


0.61988(19) 


0.79730(14) 


6.038 


0.34838(14) 


0.59472(18) 


0.76515(15) 


6.100 


0.32895(16) 


0.56187(20) 


0.72335(17) 


6.167 


0.30725(19) 


0.52423(23) 


0.67508(18) 


6.200 


0.29711(34) 


0.50662(36) 


0.65257(32) 


6.227 


0.28871(21) 


0.49280(23) 


0.63519(26) 


6.256 


0.28044(23) 


0.47834(26) 


0.61653(26) 


6.285 


0.27332(25) 


0.46636(33) 


0.60156(33) 


6.313 


0.26543(26) 


0.45313(34) 


0.58415(31) 


6.341 


0.25970(27) 


0.44252(31) 


0.57037(30) 


6.369 


0.25313(32) 


0.43176(43) 


0.55628(35) 


6.396 


0.24563(32) 


0.41907(35) 


0.54098(31) 


6.450 


0.23790(30) 


0.40540(50) 


0.52240(40) 



Table XVII; The pseudoscalar meson masses for the HISQ/tree action along the m; = 0.2?7is LCP. 
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/3 








5.900 


0.20162(09) 


0.63407(17) 


0.86972(11) 


6.000 


0.18381(37) 


0.57532(51) 


0.79046(27) 


6.195 


0.15143(14) 


0.47596(16) 


0.65506(11) 


6.285 


0.13823(50) 


0.43501(47) 


0.59951(28) 


6.354 


0.12923(15) 


0.40628(20) 


0.55982(17) 


6.423 


0.12022(12) 


0.37829(19) 


0.52161(17) 


6.460 


0.11528(21) 


0.36272(34) 


0.50137(32) 


6.488 


0.11245(15) 


0.35313(27) 


0.48716(17) 


6.515 


0.10975(12) 


0.34453(29) 


0.47516(29) 


6.550 


0.10629(16) 


0.33322(38) 


0.45989(24) 


6.575 


0.10469(68) 


0.32521(55) 


0.44869(50) 


6.608 


0.10001(17) 


0.31333(28) 


0.43286(29) 


6.664 


0.09572(18) 


0.29837(37) 


0.41178(32) 


6.800 


0.0849(18) 


0.26387(99) 


0.36257(68) 


7.280 


0.05399(43) 


0.17128(44) 


0.23766(30) 



Table XVIII: The pseudoscalar meson masses for the HISQ/tree action along the ?n; — 0.05ms LCP. 









afss 


# sources 


6.000 


0.11243(21) 


0.13224(31) 


0.15290(22) 


1 


6.195 


0.09179(21) 


0.10835(13) 


0.12525(13) 


1 


6.285 


0.08366(22) 


0.09826(13) 


0.11390(13) 


1 


6.354 


0.07825(40) 


0.09146(19) 


0.10598(11) 


1 


6.423 


0.07241(18) 


0.08515(11) 


0.09854(07) 


2 


6.460 


0.06885(11) 


0.08185(09) 


0.09454(08) 


4 


6.515 


0.06534(18) 


0.07707(15) 


0.08946(09) 


4 


6.575 


0.06104(49) 


0.07265(19) 


0.08405(14) 


2 


6.800 


0.04883(83) 


0.05774(18) 


0.06717(14) 


1 



Table XIX: Estimates of decay constants of the pseudoscalar mesons in lattice units for the HISQ/tree action along the 
mi — 0.05ms LCP. We use the normalization in which /,r ~ 90 MeV. In the last column, we list the number of source points 
used on each configuration to increase the statistics. 



p 


aMp 








6.000 


1.0292(199) 


1.1263(82) 


1.2363(24) 


1.3557(100) 


6.038 


0.9763(67) 


1.0867(40) 


1.1907(22) 


1.3277(77) 


6.100 


0.9160(117) 


1.0230(34) 


1.1211(38) 


1.2418(116) 


6.167 


0.8499(52) 


0.9570(41) 


1.0491(18) 


1.1560(55) 


6.200 


0.8300(83) 


0.9288(48) 


1.0133(24) 


1.1273(73) 


6.227 


0.8002(31) 


0.9045(37) 


0.9907(25) 


1.0662(92) 


6.256 


0.7822(27) 


0.8728(48) 


0.9548(14) 


1.0563(52) 


6.285 


0.7700(36) 


0.8602(33) 


0.9386(22) 


1.0311(25) 


6.313 


0.7462(40) 


0.8344(35) 


0.9089(21) 


1.0050(29) 


6.341 


0.7170(24) 


0.8019(30) 


0.8780(17) 


0.9705(36) 


6.369 


0.7006(21) 


0.7818(23) 


0.8572(26) 


0.9524(25) 


6.396 


0.6843(32) 


0.7605(20) 


0.8311(29) 


0.9246(37) 


6.450 


0.6505(55) 


0.7206(43) 


0.7962(18) 


0.8637(106) 



Table XX: The vector meson and the nucleon masses for the HISQ/tree action along the mi = 0.2ms LCP. 
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p 


aMp 


aMii' 




aMN 


aMn 


6.000 










1.994(13) 


6.195 


0.7562(36) 


0.8842(18) 


1.00500(93) 


1.0114(56) 


1.665(17) 


6.285 










1.4873(67) 


6.354 


0.6375(35) 


0.7499(26) 


0.85234(77) 


0.8315(95) 


1.374(60) 


6.423 


0.6047(43) 


0.6950(22) 


0.79246(83) 


0.7899(58) 




6.460 


0.578(25) 


0.6709(43) 


0.7644(22) 






6.488 


0.5647(24) 


0.6478(22) 


0.73630(68) 


0.7452(51) 




6.515 


0.5452(59) 


0.6358(40) 


0.7167(10) 


0.7247(55) 




6.550 


0.5324(24) 


0.6118(20) 


0.6929(14) 


0.7010(90) 




6.608 


0.5072(39) 


0.57572(84) 


0.65227(98) 


0.6442(69) 




6.664 


0.4732(43) 


0.5501(26) 


0.6180(10) 


0.6048(111) 




6.800 










0.8725(42) 



Table XXI: Estimates of the vector mesons, the nucleon and the Q. baryon masses in lattice units for the HISQ action along 
the mi = 0.05ms LCP. 



The integrated autocorrelation time is then given by 



T^nt « 1 + 2 ^ exp(-Ai/Ti) = 1 + 2 Y 



Ai=l 



exp(-l/Ti) 
- exp(-l/ri) 



1 + 2ti . 



(D5) 



Our data indicate that Eq. (jD3[) is a reasonable approximation for the plaquette but not for the chiral condensate 
and the topological susceptibility. To get a rough estimate of Tint we consider the data for the plaquette given in Table 
IXXIII Estimates of ti vary between 3.4-5.1 and, consequently, for Tint between 7-10. Estimates of Tint for the chiral 
condensate and the topological susceptibility could be larger due to multiple exponentials contributing long tails even 
though the Ca* are smaller. Based on these rough estimates for the plaquette, we typically save every tenth lattice 
for further measurements, for example, the quark number susceptibility. We also check for the statistical significance 
of the data for a given observable by binning to the extent justified by the statistics. 

The topological charge Q is calculated on zero-temperature ensembles following the prescription in Ref. (s^]. The 
corresponding time histories are shown in Figs. In these figures, the number of TUs examined is set to 1200 

(corresponding to the shortest time series at /3 = 7.280) to make comparison easier. 



p 


At 


□ 




Q 


6.000 


5 


0.378(42) 


0.073(14) 


0.361(57) 


6.285 


5 


0.369(29) 


0.076(19) 


0.262(39) 


6.460 


5 


0.336(28) 


0.120(21) 


0.397(51) 


7.280 


6 


0.168(34) 


0.025(40) 





Table XXII: Estimates of the autocorrelation coefficient Ca* with At in TU for a representative set of the HISQ/tree zero- 
temperature lattices measured using the plaquette, strange quark chiral condensate and the topological charge. 
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Nr=9, 


iVr = 12 


p 


□ 






□ 






□ 






6.000 


0.172(18) 


0.119(14) 


0.219(22) 


- 


- 




- 


- 


- 


6.100 


0.158(09) 


0.219(06) 


0.270(11) 


- 


- 


- 


- 


- 


- 


6.195 


0.121(16) 


0.313(29) 


0.371(23) 


0.156(19) 


0.096(18) 


0.226(29) 


- 


- 


- 


6.285 


- 


- 


- 


0.116(09) 


0.163(12) 


0.218(15) 


- 


- 


- 


6.354 


0.085(12) 


0.068(18) 


0.295(39) 


0.168(34) 


0.262(33) 


0.328(29) 


- 


- 


- 


6.460 


- 


- 


- 


0.155(09) 


0.353(21) 


0.413(30) 


- 


- 


- 


6.550 


0.129(16) 


N/A 


0.172(21) 


0.108(17) 


0.238(37) 


0.413(36) 


- 


- 


- 


6.664 








0.132(29) 


0.078(20) 


0.337(32) 


0.117(37) 


0.285(34) 


0.342(42) 


6.800 


0.027(17) 


N/A 


N/A 


0.117(13) 


0.043(10) 


0.209(10) 


0.091(22) 


0.392(46) 


0.298(42) 


6.950 








0.063(28) 


N/A 


0.085(14) 


0.121(12) 


0.355(31) 


0.385(41) 


7.030 














0.116(21) 


0.278(32) 


0.292(34) 


7.150 














0.091(04) 


0.132(29) 


0.177(15) 



Table XXIII: Estimates of the autocorrelation coefficient Ca* with At = 10 for a representative set of HISQ/tree finite 
temperature lattices measured using the plaquette, and the light and strange quark chiral condensates. The correlations 
decrease away from the transition region, therefore in some cases, we were unable to extract the autocorrelation coefficient. 
Note that the chiral condensate is calculated with stochastic estimators that introduce additional noise in the correlation 
function. 
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Figure 30: Time history of the topological charge measured every 5 time units on 32* configurations at /? = 6.285 with the 
HISQ/tree action. 
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Figure 31: Time history of the topological charge measured every 5 time units on 32^ x 64 configurations on /3 = 6.460 with 
the HISQ/tree action. 
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Figure 32: Time history of the topological charge measured every 6 time units on 48^ x 64 configurations at 13 — 7.280 with 
the HISQ/tree action. 



[46] A. Bazavov et al. (MILC collaboration), Phys. Rev. D82, 074501 (2010), 1004.0342. 
[47] S. R. Sharpe, PoS LAT2006, 022 (2006), hep-lat/0610094. 
[48] M. Creutz, PoS LAT2007, 007 (2007), 0708.1295. 

[49] U. M. Heller, F. Karsch, and B. Sturm, Phys. Rev. D60, 114502 (1999), hep-lat/9901010. 

[50] P. Hegde, F. Karsch, E. Laermann, and S. Shcheredin, Eur. Phys. J. C55, 423 (2008), 0801.4883. 

[51] M. Clark, A. Kennedy, and Z. Sroczynski, Nucl. Phys. Proc. Suppl. 140, 835 (2005), hep-lat/0409133. 

[52] M. Clark, P. de Forcrand, and A. Kennedy, PoS LAT2005, 115 (2006), hep-lat/0510004. 

[53] A. Bazavov et al. (MILC Collaboration), PoS LATTICE2008, 033 (2008), 0903.0874. 

[54] R. Sommer, Nucl. Phys. B411, 839 (1994), hep-lat/9310022. 

[55] C. Aubin, C. Bernard, C. DeTar, J. Osborn, S. Gottlieb, et al., Phys. Rev. D70, 094505 (2004), hep-lat/0402030. 
[56] S. Necco and R. Sommer, Nucl. Phys. B622, 328 (2002), hep-lat/0108008. 
[57] D. Toussaint, private communications. 

[58] A. Bazavov et al. (MILC Collaboration), PoS LATTICE2010, 074 (2010), 1012.0868. 

[59] C. Davies, E. FoUana, I. Kendall, G. Lepage, and C. McNeile (HPQCD Collaboration), Phys. Rev. D81, 034506 (2010), 
0910.1229. 

[60] A. Gray, I. Allison, C. Davies, E. Dalgic, G. Lepage, et al., Phys. Rev. D72, 094507 (2005), hep-lat/0507013. 
[61] C. W. Bernard, T. Burch, K. Orginos, D. Toussaint, T. A. DeGrand, et al., Phys. Rev. D64, 054506 (2001), hep- 
lat/0104002. 

[62] W.-J. Lee and S. R. Sharpe, Phys. Rev. D60, 114503 (1999), hep-lat/9905023. 

[63] M. Cheng, N. Christ, C. Jung, F. Karsch, R. Mawhinney, et al., Eur. Phys. J. C51, 875 (2007), hep-lat/0612030. 

[64] S. Ejiri, F. Karsch, E. Laermann, C. Miao, S. Mukherjee, et al., Phys. Rev. D80, 094505 (2009), 0909.5122. 

[65] J. Engels, S. Holtmann, T. Mendes, and T. Schulze, Phys. Lett. B492, 219 (2000), hep-lat/0006023. 

[66] D. Toussaint, Phys. Rev. D55, 362 (1997), hep-lat/9607084. 

[67] J. Engels and T. Mendes, Nucl. Phys. B572, 289 (2000), hep-lat/9911028. 

[68] J. Engels, S. Holtmann, T. Mendes, and T. Schulze, Phys. Lett. B514, 299 (2001), hep-lat/0105028. 

[69] O. Kaczmarek, F. Karsch, E. Laermann, C. Miao, S. Mukherjee, et al., Phys. Rev. D83, 014504 (2011), 1011.3130. 

[70] A. Cucchieri, J. Engels, S. Holtmann, T. Mendes, and T. Schulze, J. Phys. A35, 6517 (2002), cond-mat/0202017. 

[71] D. Wallace and R. Zia, Phys. Rev. B12, 5340 (1975). 

[72] P. Hasenfratz and H. Leutwyler, Nucl. Phys. B343, 241 (1990). 

[73] A. V. Smilga and J. Stern, Phys. Lett. B318, 531 (1993). 

[74] A. V. Smilga and J. Verbaarschot, Phys. Rev. D54, 1087 (1996), hep-ph/9511471. 

[75] M. Cheng, S. Datta, A. Francis, J. van der Heide, C. Jung, et al., Eur.Phys.J. C71, 1564 (2011), 1010.1216. 
[76] RBC-Bielefeld Collaboration work in progress. 

[77] O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow, Phys. Lett. B543, 41 (2002), hep-lat/0207002. 
[78] S. Digal, S. Fortunato, and P. Petreczky, Phys. Rev. D68, 034008 (2003), hep-lat/0304017. 
[79] O. Kaczmarek and F. Zantow, Phys. Rev. D71, 114510 (2005), hep-lat/0503017. 
[80] P. Petreczky and K. Petrov, Phys. Rev. D70, 054503 (2004), hep-lat/0405009. 

[81] C. Bernard, C. E. DeTar, L. Levkova, S. Gottlieb, U. Heller, et al., PoS LAT2007, 090 (2007), 0710.1118. 
[82] A. Bazavov et al. (MILC collaboration), Phys. Rev. D81, 114501 (2010), 1003.5695. 



